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Department of the Inteutou, 
% IT. S. Geolooioal Survky, 

Washington, 7>. 0., May lily 1887. 

Sir: I have the honor to trausinit herewith the results of certahi 
investigations, which may be broadly designated as relating to the Ibnu 
and position of the sea level. These investigations wore begun in part 
previous to my connection with the Geological Survey, but they wore 
taken up again in 1885, with your approval, at the request of Mr. (h K. 
Gilbert and Prof. T. C. Chamberlin, for solutions of some special prob- 
lems which arose in their geological researches. Tiie work has boon 
prosecuted simultaneously with other lines of ofUce and field work. It 
reached its present form substantially, however, inore than a year agoj 
and the principal numerical results of the discussion of ProfesHor Cham- 
berlin's problem are incorporated with his paper 0!i The DriftloMH 
Area, in the Sixth Annual Report. The purely mathematical featuroH 
of the paper have been published also in the Annals of MathcmaticH^ 
Nos. 5 and 6, vol. 2, and No. 1, vol. 3. I have delayed ottering the 
complete manuscript for publication up to this time in onh*.r that I 
might give it a careful revision and check all the more important form- 
ulas by independent processes of derivation. 

The qnestions treated in this paper are for the most part tunutHHiivily 
somewhat mathematical. They are, however, fundamental quimtiotm in 
geophysics, and although the mathematical form of presentation han 
been followed throughout, an attem]it has been made to sUite the end 
results and formulas in such a way that they may l>e understood and 
used with safety by those who may not care U> follow the detaiU of the 
analyslA. For the benefit of such readers a key U) the matbematii^al 
symbols employed is given in addition to the liMt of cA}nUtntH and khU' 
eral index. 

WhOe the analysis of this investigation was design^e^l €tH\H*ii\ii\\y to 
solve the particakur problems of 3fessrs« Gill^^f t and Chamfierlin, ft hm 
not been confined to those probleuM, bat ban \Htitn a^lapte^l to ihn itni\nt 
class of problems to which they belong. It hi UoiHtA^ Iherefore, that 
the resnlts of the paper will be of interest ami value to ffiffMlnnlHUt and 
mathematicians as well as to geoloj^ts. 

Very respeetfallj, yoar obedient servant, 

Hon. J. W. Powell, 

IHreet^r U. A Geolofhal Surrey. 



ON THE FORM AND POSITION OF TEE SEA LEVEL 



By E. 8. Woodward. 



I. INTRODUCTION. 

1. The problem of the form and dimensions of the sea level surface 
of the earth has been one of peculiar difficulty. The combined efforts 
of the ablest mathematicians of the past two centuries, supplemented 
by the most laborious and costly geodetic measurements have yielded 
us the first approximation only to the complete solution. Fortunately 
this first approximation is exceedingly close. It assigns to the sea level 
a form which differs but slightly from that of an oblate spheroid, whose 
major and minor semi-axes are about 20,926,000 and 20,855,000 English 
feet, respectively. This spheroid, or reference ellipsoid, as it is some- 
times called, has its minor axis coincident with the earth's axis of rota- 
tion and is usually regarded as sensibly fixed in position and dimen- 
sions. With respect to it the actual sea surface or geoid must be 
imagined to lie partly above and partly below by small but unknown 
amounts, the determination of which, if possible, will constitute a 
second approximation to the figure of the earth. For many if not most 
of the applications of science the reference ellipsoid suffices^ the first 
approximation is nearly enough correct. But geodesy, on the one hand, 
has attained such a degree of perfection in precise measurement that 
the discrepancies now brought to light in some of its operations must be 
attributed largely if not chiefly to defects in theory. These discrepan- 
cies must be explained before any considerable advance can be expected 
in our knowledge of the figure of the earth along the present lines of in- 
vestigation. Their true explanation is apparently intimately connected 
with the form of the geoid, and it is to the study of the form, therefore, 
rather than to the determination of the dimensions of the geoid that we 
may look for future progress in geodesy. Geology, on the other hand, 
has raised many questions relative not only to the form, position, and 
fixity of the geoid proper, but also with respect to the allied equipoten- 
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16 FORM AND POSITION OF THE SEA LEVEL. [bull. 48. 

tial surfaces of isolated bodies of water at higher or lower levels. It is 
found, for example, in geological investigations, that the shore lines of 
extinct seas do not always coincide with existing level lines, but often 
cross them at decided angles, or that the water level lines traced on 
islands in such extinct seas differ in elevation from contemporaneous 
lines traced on their distant shores. Aside from the changes which 
may have been due in these cases to subsidence or upheaval, the ques- 
tion may be raised whether such slopes or differences in elevation rela- 
tive to present level lines may not have been caused by adjacent 
attracting masses, which have since disappeared, like the ice mass of 
the glacial epoch, or in a lake basin by the presence of the water itself. 
Correct and complete answers to such questions require a knowledge 
of the existing geoid and of the causes which may have produced secu- 
lar variations in its form and position. 

At present it is by no means clear how any specially extensive, addi- 
tions to our information concerning the more minute features of the sea 
surface are to be obtained. It may even be doubted whether we have 
not reached a practical limit in the first approximation ; whether in fact 
the distribution of matter within the earth's surface is not so irregular 
as to preclude gaining anything more than an empirical formula for the 
deviations of the geoid from the ellipsoid of reference. It seems prob- 
able, however, that the forces producing these deviations have their 
seat in a comparatively thin terrestrial crust resting on a fluid or plastic 
substratum (or nucleus), or that such was the antecedent condition of 
the earth, and that our failure to perceive the relations of the crust to 
the substratum is the chief obstacle to improvement. 

In the absence of a complete rational theory the best evidence which 
analysis can bring to bear on questions pertaining to the geoid is largely 
of a negative character. The effects which would result under certain 
conditions can be computed, but it is not always possible to prove that 
those conditions accord with the actual facts. Investigation must 
proceed to some extent upon doubtful postulates, and computations 
must be made from uncertain data. But notwithstanding this limita- 
tion -on the calculations we are about to consider, they will generally 
possess a value in excluding or confirming hypotheses, or in furnishing 
limiting values for the eflfects of observed causes. 

2. A considerable class of problems concerning the sea level is that 
in which the attracting or disturbing mass is symmetrically disposed 
about a radius of the earth's surface, and is situated on or near the 
surface. As examples of this class we may adduce the two following, 
which led to this investigation : 

(a) Given the dimensions of a lake basin having a circular border. 
When the lake was full of water it left a trace of its surface along the 
border and on an island at the center of the basin. After the water 
had disappeared a line of spirit levels was run between the water tr^roe 

(100) 



wooDWABD.] b.esum:6 of results. 17 

on the island and that on the border ; what difference in altitude shoald 
have been fonnd V 

{b) Assuming the accumulation of ice in glacial times to have been 
in the shape of a spherical stratum bounded by a circle, or some sort of 
meniscus symmetrical about an axis, and that the earth's crust did not 
yield under the weight of the ice, what were the resulting distortions 
in. the sea level P 

It will be seen that these problems are essentially the same. They are 
substantially identical also with the problem of the effect of continental 
masses on the sea level, since the continents may be represented, ap- 
proximately at least, as spherical strata having circular borders, or as 
masses of meniscoid shape. 

3. The following paper is devoted to the investigation and discussion 
of this class of problems. An attempt has been made to develop the 
theory of their solution so far as is necessary to render practicable the 
numerical evaluation of the characteristic effects of the disturbing mass 
in any special case. In Articles II to XI the theory of the effect of a 
mass in the shape of a spherical stratum having a circular border and 
uniform thickness is worked out with considerable detail. The only re- 
strictions imposed on this mass are that its density is uniform, and that 
the ratio of its thickness to the earth'js radius may be neglected in com- 
parison with unity. Expressions for the potential of the disturbing mass 
at any point of the disturbed surface are derived in terms of a definite 
integral and in terms of spherical harmonics; and the degree of approx- 
imation of these expressions is investigated. Equations to the dis- 
turbed surface are assigned for the case in which the effect of the re- 
arranged free water is considered, as well as for the case in which that 
effect is neglected. The disturbance in the former case is shown to be 
equal to that in the latter, which is expressed in compact integral form, 
plus a rapidly converging series of additive terms. 

' This problem was proposed by my colleague, Mr. G. K. Gilbert, to tlie mathematical 
section of tlie Washington Philosophical Society, February, 1884. In his geological 
investigations within the area of the Qnaternary sea known as Lake Bonneville, Mr. 
Gilbert has found traces of the central portions of the ancient lake surface to be more 
than 100 feet higher than the traces of the contemporaneous surface at its margin. 
A complete consideration of the effects of the causes which might contribute to this 
distortion requires, obviously, a numerical evaluation of the depression of the level 
surfaces within the area, due to the removal of the water. 

A solution of the problem was given by the writer before the above-named society 
in March, 1884, and a more complete discussion will be found in sections 39-41. 

2 To what extent the form and position of the sea level may be modified by the mere 
attraction of glacial masses is a question which ha.s been much discussed by geologists, 
it was proposed to the writer by Prof. T. C. Chamberlin, geologist in charge of the 
division of glacial geology, U. S. Geological Survey. The question is considered at 
some length in sections 42>52, and a review of the work of the more prominpnt 
mathematicians who have discussed the problem is given in sections 53-62. The 
principal numerical results of the writer's investigations are given in Professor 
Chamberlin's paper on The Driftless Area, in the Sixth Annual Report of the U. S. 
Geological Survey. 
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18 FORM AND POSITION OF THE SEA LEVEL. [bull. 48. 

In Article XII the Investigation is extended so as to assign the eflfect 
of any mass of uniform density having a symmetrical distribution aboat 
a radius of the earth's surface. Particular attention is paid to a class 
of masses whose shapes are assigned by a formula which represents 
fairly well the mass features of the problems (a) aud (b) above. 

Under the head of applications, Articles XllI to XVI, the character, 
istic properties of the equipoteutial surfaces in a lake basin are first 
considered. Then the variations in sea level attributable to continental 
glaciers or ice caps are discussed at some length. The angular radial 
extent of the ice mass is, for the most of the discussion, assumed to be 
380, for the reason that this is the extent of a mass of nearly uniform 
thickness, which would produce the maximum upheaval of the water 
along its border. The external shapes of the various masses, their vol- 
umes, and the distortions of the sea surface attributable to them are 
given in detail. The minimum thicknesses of ice masses of varying 
angular radial extent, requisite to produce av^erage slopes of five feet 
per mile within one degree (69 miles) of their borders, and the extent of 
variation in sea level on the hypothesis of an alternation of glaciation 
at the poles of the earth, are also worked out. 

In the historical note of Article XV, the allied investigations of Arch- 
deacon Pratt, Mr. D. D. Heath, and Sir William Thomson, on the prob- 
lem of glacial submergence are reviewed. The special cases they have 
considered are shown to be easily derived from the general formula of 
Article XII. 

Finally, in Article XVI, a brief discussion of the eflfect of continental 
masses in distorting the sea level is given. It is shown that according 
as the continents are or are not superficial masses unbalanced in their 
attractive effects, the sea surface must be very irregular or deviate only 
by minute quantities from the ellipsoidal form. It is also shown that 
although a continent whose radial element masses are in a condition 
bordering on hydrostatic equilibrium would produce but slight disturb- 
ances in the position of the sea level, it might nevertheless cause a.con- 
siderable slope of the sea surface, or.deflection of the plumb line along 
its border. 

A. THEORY. 
II. MATHEMATICAL STATEMENT OF PROBLEM. 

4. The solution of the general problem outlined in the preceding sec- 
tion depends on the principle of hydrostatics that the potential of the 
forces producing a liquid surface in equilibrium has a constant value 
for all points of that surface. In the case of the earth, if the potential 
of all the attractive forces acting on a unit mass at any point of the sea 
surface be denoted by P, the distance of the point from the earth's axis 
of rotation by 15, and the velocity of rotation by a?, the form of the sar- 
face will be completely defined by the equation 
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P+iPcD^=zQ, constant. (1) 

The exact value of P in this equation is a complicated function of the 
densities of the element particles of the earth and of the co-ordinates of 
those particles and the attracted point. For the present purposes, how- 
ever, it will be sufficient to consider P due to a centrobaric sphere of 
equal mass and volume with the earth and concentric with the earth's 
center of gravity. Since we shall only consider relative positions of 
any point on the sea surface, the potential due to centrifugal force, which 
is represented by the second term in (1), may be neglected. 
^ 5. If Oo and 6© denote the equatorial and polar semi-axes, respectively, 
of ther earth's ellipsoid, and Tq the radius of the sphere just referred to, 

n= 'y/~ao%. X2) 

Using Clarke's values^ of ao and bo we have 

ao=20926062 English feet, 
6o=20855121 English feet, 
ro=20902394 English feet, 
log ro=7.32020. 

The surface of the sphere thus defined may be regarded as the sur- 
face assumed by a thin film of sea water covering a nucleus whose mass, 
plus the mass of the film, equals the earth's mass. We shall call this 
ideal surface the undisturbed surface. With respect to it the real sur- 
face of the earth lies partly without and partly within ; but so far as 
small relative changes in sea level are concerned it is practically imma- 
terial whether we refer to the actual closely spheroidal surface or to 
the simpler spherical one. 
6. Let 

ilf= mass of the earth, 
/)„=mean density of earth. 
Then, 

M=^7tro^p^, (3) 

and the equation to the undisturbed surface is 

^=^;rroV.= Oi, (4) 

Vo o 

Ci being a constant. 

Suppose, now, a new mass, «i, of density p (positive or negative) be 
placed in any fixed position relatively to the undisturbed surface. The 
resulting sea surface will then differ from that defined by (4). To de- 
termine this difference let V be the potential of the disturbing mass m 

* Comparisons of Standards of Length, made at the Ordnance Survey Office, Soath- 
ampton, England, by Capt. A.. R. Clarke, R. E. Pablished by order of the secretary 
of state for war, 1866. 
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at any point of the disturbed surface, and let v denote the elevation or 
depression of this point relative to the undisturbed surface. The equa- 
tion to the disturbed surface will then be 

M 
—— + V= O2 , a constant. (5) 

n+v 

The difference of this and (4) to terms of ihe first order inclusive in 
V is 

whence, putting 

v=(V^Vo)^. (6) 

Since M/rQ^=g, the velocity increment at the earth's surface due to 
the earth's attraction, (6) may be written 

v=l — -\ (6y 

9 

Vo in the last two equations is the value of V when 17=0, or the value 
of F along the line of intersection of the disturbed and ^undisturbed 
surfaces. If we put 

-+-o=y|=I (7) 

This equation represents the elevation of the disturbed surface above 
a spherical surface of equal potential, whose value is 

M _^ 



n—vo 



since the difference between this and (5) gives (7). 

The constant Vo may be determined from the obvious condition that 
the disturbed and undisturbed surfaces must contain equ<al volumes. 

It is evident that the equations just derived will hold true if the 
mass m be a part of the earth's mass, so long as the ratio m/M may be 
neglected relatively to unity. Thus, in the problems we shall consider, 
m may represent the mass of a continent, the deficiency in mass of a 
lake or lake basin, or the ice mass of the glacial epoch. 
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III. EVALUATION OF POTENTIAL F— DISTURBINa MASS ON SURFACE 
OF EARTH, OF UNIFORM THICKNESS AND DENSITY, AND WITH 
CIRCULAR BORDER. 

7. The next step in the solution requires the determination of the 
potential Fof the attracting mass for any point of the disturbed sur- 
face, whether without or within the circle which we have assumed to 
define the boundary of the mass. Although the nature of the mass may 
be such as to prevent the water from permeating it freely, the surface 
the water would take if not so restricted is an essential part of the 
disturbed surface. 

la order to derive an expression for F, let the rectangular and polar 
co-ordinates of any point of the attracting mass be defined by the usual 
relations, viz : 

x=r cos 6 cos A, 

y=:r cos sin A, 

z=zr sin 6, 

in which and A correspond to polar distance and longitude, respect- 
ively, the position of the origin being arbitrary. With reference to 
the same origin, let the co ordinates of the attracted point on the sea 
surface be 

ir'=r' cos 6^ cos A', 

yf=r^ cos 6' sin A', 
2/=r' sin 6\ 

If D denote the distance between the attracting and attracted points 
and 

cos ^=cos COS ^'+sin 6 sin 6^ cos (A— A'), (8) 

2)2=r2+r'2— 2rr' cos ^=(r-r')^+4rr' sin^ t . (9) 

The volume element of the attracting mass is 

dxdydz=:Y^dr sin OdOdX, 

Hence, if p denote the density of the attracting mass, a general ex- 
pression for the required potential is 

n p rrHr Bin dd6d\ ,-^. 

^='J J J — 5 ^''^ 

We must now evaluate this integral. Taking the center of the sphere 
of reference as the origin of co ordinates, let 

r=ro+w, 
and (11) 
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in which u and v are small quantities relative to ro, v being the same 
of course as defined by equation (6). Premising what will be proved in 

Article IV, namely, that we may neglect quantities of the order — , — , 

2 



and (^ — ^ ) equation (9) gives 



i)=2ro sin ^. (12) 

Prom the first of equations (11) 

and (13) 

to terms of the order - . 

ro 

As to the magnitude of the quantities neglected, it may be remarked 
in passing that ro is in round numbers 21,000,000 feet (see section 5), 
while u and v may be restricted to values less than 100,000 feet; so 
that the fractions neglected will not exceed s^^ . 

Without loss of generality we may assume the line from which and 
0^ are reckoned to pass through the attracted point, and the plane from 
which A and X' are reckoned to pass through the attracted point and 
the center of the attracting mass. In this case ^'=0 and X'=0, and (8) 
gives tp=0. 

By means of this relation and the equivalents in (12) and (13) the in- 
tegral in (10) becomes 



=^"^/// 



du cos ~ dddX. (14) 



If the uniform thickness of the attracting mass be denoted by A, the 
limits of w in (14) will be and h. Let the limits of 0, which are ob- 
viously functions of A, be denoted by 0i and 02. The limits of A are 
evidently equal in magnitude but of oppdsite signs. Hence we have 

V=2roP rdu(\o^^^d0 rdX^^nhp JY Bin J- sin ^^d\. (15) 

8. To complete the evaluation of (15) it will be convenient to change 
variables. Consider the spherical triangles formed by the attracted and 
attracting points, the center of the attracting mass, and the points in 
which the arc cuts the circle bounding the mass. Thus in Figs. 1 and 
2, let P be the attracting and A the attracted points, C the center of 
the attracting mass, and BDO the bounding circle. Then 

0=AP2LnaX=BAC. 
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Draw OE perpendicular to AB and put 

CE=p, 



From either figure 



6=q+8, 



6^]=^— «o, 



whence 



sm ^— sin^=2 



BG=/3, 
BE=8o, 
AE=q, 



cos ? sin ^ . 



^2=g+«o, 



(16) 



The right-angled spherical triangles of either figure give 



«^« ^ cos a ^,^„ 

cos (7= , cos 

cos p 



cos S • • •1 

«o= ) sin ©=sin a sm A. 

cos 2) 



(17) 





Fio. 1. 



Fig. 2. 



The first two of (17) give 



2cosS^=l+^«-^, 
2 cos p 



whence 



2sin2?9=l-2£?^, 
2 cos^' 



2 cos ^ sin ^o^f(cosjp+c os Qr)(cosjp-co s/?)]^ 
2 2 cos^ 



(18) 



From the last of (17) 



(cos^l?— cos^a)*' 



(19) 



Now, the last of equations (17) and the diagrams show that the limits 
of jp, corresponding to the limits of A, are and or or and /?, according 
as the attracted point is within or without the circle bounding the at- 
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tracting mass. Hence, if we denote the potentials in the two cases by 
Vi and V2, respectively, the equivalents in (15), (16), (18), and (19) give 



i / cos J?— cos /? \^, 
j Vcosj?— cosay ^' 



T^'=4roAp -"_Z" 4p, (20) 



«5/?; 



F.= 



=4roftpfV?2?^=52i^Y^i'. ^ f21) 

^ I V^osj)— cos aj ^^ ' ' 

9. The integrals in these equations are in general elliptics of the third 
species. They may be evaluated by the usual processes applicable to 
elliptics, by series, or by mechanical quadrature. 

The integral in (20) presents some apparent difficulty, since the ele- 
ment function is infinite at the upper limit, except when «=/?. Again, 
in case a=0, this integral assumes the anomalous form 



j](i=3^'* 



6 
the value of which is n sin tL as may be easily verified by means of (15), 

(16), and (18). These peculiar features may be removed by the follow- 
ing change of variables, which secures the same constant limits for both 
(20) and (21). 



For brevity put 




J J /C08ff-C08/?Y^ (22) 

■ \cos2)— cosay '^^ ^ ' 



Then, observing that 



put in Ji 



and in Z 



2 



^ I / C0Sj?-C08/g Y^p 

' j \cos2>— cosay 



. sin^^-sin^S 
cos 2>— cos /?__ 2 2 

cos »— cos a"~ . 9 a • o p ' 
■^ siu^ ^— sm^^ 



sin^=sin-sin;^i, 

sin:^=sm^ sm y^^ 
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These give 

dp= J 

(l-sin^^sin^;/,)* 



25 



and 



dp= 



2 sin ^ cos y2dy2 



(I— sin2^sin^7^2)* 
and the limits for both yi and y2 are and ^. Therefore 



2 



2 



2i= 



2sin|i 



sm* — 



1 ^ ■ 2 



sin-1 




dyi 







(1— siii^^sin- ;^i)* 
a</3',. 



(24) 



2 



r.= 



2 sin^ -^ cos^ ^2^X2 




.2^ 



(25) 



i , 



-_ |sinV2 I ( l-sin2,^sinV2^ 



sin-- 



ft 

10. Some special values ot the integrals (22) to (25) and the corre- 
sponding potentials (20) and (21) are worth deriving. These values are : 

(a) For a point of the disturbed surface at the center of the disturb- 
ing mass, a=0, and (24) gives 



/3 



2Jo 2 



(26) 



and the corresponding value of the potential is 

Fi =4ro/^p7r sin P-. 

(h) For a point at the border of the disturbing mass, a 
from (24) and (25) 



(27) 
■/3, and hence 



Ii=Ii=2 arc sin Tsin ^ sin ^ j '=A 

(109) 



(28) 



y=0 
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a result which is reached more readily from (22) or (23). The corre- 
spondiDg value of the potential is 

Vi=V^=^rohp/3. (29) 

(c) For a point of the disturbed surface 180° from the center of the 
disturbing mass, a=7rj and (25) gives 




IT 

Sin2 1 C0S2 y2dy2 /^^(l -COS^I^d;/^ 



/.=2 / f _^ 



1 — sin2 P gin2 y^ I ^j^ QQ^2 P ^^ y^ 

Z 9/0 2 



=;r— 2co8 ^ arc tan 



rtan( cos^tan;/2 j 1^ 



= Trfl - cos|^ =2;r sin^ |, (30) 

and the potential is 

V2=Srohp7r&m^^. (31) 

{d) Suppose y^=7r ; in other words, let the attracting mass cover the 
whole sphere. Then (24) gives 

Ii=7r, (32) 

and^e have the following well-known approximate value for the po- 
tential of a spherical shell for a point on its interior or exterior surface^^ 

viz : 

^VohpTT. (33) 

This result follows also from (27) or (29) if we make /J=;r. 

IV. DEGREE OF APPROXIMATION OF THE EXPRESSIONS FOR THE 

POTENTIAL V. 

11. In deriving the expressions (20) and (21) for the potential of the 
disturbing mass, it was assumed that a sufficient degree of approxima- 
tion is attained if quantities of the orders - , - , and upwards are 

To To 

neglected. The grounds of this assumption need to be examined with 
some care. For this purpose we shall derive the exact expressions, in 
form at least, for the potential of the disturbing mass at its center, at 
its border, and at 180^ from its center. A comparison of these exact 
values with the approximate values given by (27), (29), and (31) will 
show the order o? approximation of (20) and (21). 

1 See equation 34. 

(110) 



WOODWARD.] DEGREE OP APPROXIMATION ATTAINED. 27 

We will first write down the expressiou for the potential of a com- 
plete spherical shell, of uniform thickness and density, for a point within 
its bounding surfaces. This expression will be useful, as a check on 
formulas relating to partial shells. 

Let the radius of the interior surface of the complete shell be r©, the 
uniform thickness A, the uniform density p, and the distance of the at- 
tracted point from the center of the shell Tq+v. 

Then the potential is^ 

F=4roAp;r(^l+^^^,-p^-2(^^ . (34) 

12. Let the notation be the same as that used heretofore ; i. e., let r 
be the radius- vector of any point of the attracting mass, r' the radius- 
vector of any point of the disturbed surface, p the density, h the thick- 
ness, and /3 the angular extent of the mass; and ff the angular distance 
between the attracted and attracting points. For points of the dis- 
turbed surface lying above the undisturbed surface, r will be less than 
r'=ro+v over the range r© to r'y and r will be greater than r' over the 
range r* to Vo+h. Bearing these facts in mind it follows that the exact 
Value of the potential of the disturbing mass for the point where its 
axis pierces the disturbed surface is 

^ ^ . o , M sin 6d8 

I I vr^+r'^— 2rr' cos 6 



, , sin 6d(^ 

^^""X '-"^'- \ y/^+r'^-2rr' cos d 




ro4-o 



l^o+w 



=2p;r I ^^4rr'sin2|+(r'-r)2-(r'-r) Y^J-+ 



fcrn+Zt 



2p;r 1 ('-/4rr'sin2.-|+(r--r')2--(r-.rO^-^. 



See Price^s Calculus, vol. 3, p. 299. 
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Since 
j4»-r'sm»^,+(r-r'f=28in^VFp/l+-^^^^=^-- • A = 

we find by expansion, integration, and reduction, to terms of the first 
order inclusive, 

3Asin 7j —A 

1+ — ^ 



4(ro+i?)sin^ 
F= 4ro/^p;r sin^^ ^. (35) 

2(ro+^) sin '- 2{ro+v)h sin ^ 

If we make /3=7t this expression agrees with (34) to terms of the 
second order. 

Equation (35), it will b^ observed, differs from (27) by certain terms 
which must be small unless h and v are very large. In one of tlie most 
important applications discussed in the sequel, ^=10,000 feet, 2?=3,000 
feet, and /3=SS^. With these values, since u is in round numbers 
21,000,000 feet, the quantity within the parentheses of (35) differs from 
unity by less than yoioo* I* )5=60o, which is (see section 23) the an- 
gular extent of mass required to produce the maximum elevation of the 
disturbed surface at the center of the mass, the quantity within the 
brackets of (35) exceeds unity by less than 8oVo» using the above values 
of h and v. 

13. Similarly, if the attracted point be at the border of the attracting 
mass, the exact valile of the potential is 



2 



1 ^ 



+^/^ / p^*' / J(r-rOH4rr'sin2^ • ^'^ 

r« *- r.+v 
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in which 

siii2 ^= s iu^/?co s^/\ ^ 
2 i— sin^/^siu^A* 

Introducing this value, the first of the above integrals relative to \ 
becomes 




;iu2 6 cos2 A + f!l_L Vcos^ /? 



Since the numerator of the element fuQction in this integral is greater 



r'—r 



than sin /? cos A and less than sin /? cos AH — ,— — cos/?, the value of the in- 
tegral lies between (r'+r)/5 and 

IT 

Suppose the exact value of tbis integral is 

(r'+r)/j|l+(?.^|. 

Likewise, represent the exact value of the second integral relative to 
A. by 

(r+r')^|l+<;),^|. 

Then to terms of the first order inclusive the above expression for 
the potential becomes 

The first term of this agrees with (29). 
The quantities Qx and Qi lie between 



and 




cos /? / ^A 



• _ _- 



/? / \/l-sin2/^sin2A' 



When /?=n> or when the attracting stratum covers a whole hemi- 
sphere, ^1=^^=0, and (36) becomes 



^^=2r«/tp;r 1 1+..,--,— ,- .„ '^\ \u \ • (36') 

«■ f -^ ^2(ro+v) 2(ro+v)h } ^ ' 

(113) 
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This agrees, as it should, to terms of the second order with the half of 
(34). For ^=10,000 feet and i? =3,000 feet the quantity in the paren- 

thesis exceeds unity by less than :4^o7y. When /?<^, Qi and Qgwill ex- 

ceed 0, and the degree of approximation of (36) will be somewhat higher 
than that of (36'). 

14. For a point of the disturbed surface ISO^ from the center of the 
disturbing mass the exact value of the potential is 

e /ro l/o Vr2+r'2+2rr' cos 6 

{r+r'— y/~t^-\-r'^+2rr' cos yS)^^^. 



When 



(r— r')^<4rr'cos2^ , 



which is the only case we need consider, 



y/r^^r''i+2rr' cos /^=2 cos ^^^'^'^'X ^+ — ^ ^~ • • • 



= 2rcos^)l+-(^""^')' 



^ ' 8rr' cos ^ 



Srr'cos^ 



r—r' i 



Hence, expanding, integrating, and reducing, there results to terms of 
the first order inclusive, 

F=8.o/^p.8in^|{l+^^^;, (37) 

the first term of which agrees with (31). Using the values ro=21,000,000 
feet, ^=10,000 feet, and t?=— 1,000 feet (v being here intrinsically nega- 
tive), the factor in the brackets of (37) exceeds unity by 



1 



2600 • 



V. DEVELOPMENT OF POTENTIAL V IN SERIES OF SPHERICAL 

HARMONICS. 

15. The preceding expressions for the potential of the attracting mass, 
namely, F, as defined by equations (20) and (21), are sufficient for most 
of the applications to bo considered in the sequel. They possess the 
obvious advantage of a compact integral form. For some purposes, 
however, it will be desirable to have V expressed in a series of spheri- 
cal harmonics or Laplace's functions. We may thereby arrive at equa- 
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tions (20) and (21) by a process diflfering from that followed in Article III 
and establish a useful harmonic development of the elliptic integrals 
Ji and I2 of equations (22) to (25). 

16. Expressions fulfilling the present requirements may be derived 
from equation (10) by expanding D-^ in a series of ascending powers of 

-J and — . Thus from equation (9) we have 

♦ B=p[-Po+Pi(p)+a(J) + ---] wlienr<r', 

i=i[Po+Pi(0+P.(^y+ . . . ] when r>r'. 

In these equations Po, Pi, P2, etc., are Laplace's coefficients of the 
zero, first, second, etc., order, respectively. They are functions of the 
angular coordinates ^, ^, A, and V only. 

Taking, now, the center of the sphere of reference as the origin of co- 
ordinates, and supposing the line from which 6 and 6' are reckoned to 
pass through the center of the attracting mass, we shall have for a mass 
of uniform thickness /t, all integrations ia (10) independent. For that 
part of the disturbed surface which lies above the undisturbed surface 
r'^zrQ+v will fall between the extreme values of r, which are n and 
ro+A; and hence r<^r' for values between Vo and n+v^ r^r' for values 
between ro+t? and ro+/^ For that part of the disturbed surface lying 
below the undisturbed surface r>r'. In both cases the limits of 6 are 
0, and the angular radius /? of the attracting mass and the limits of A 
are and 2;r. Therefore for that part of the disturbed surface lying 
above the undisturbed surface equation (10) gives 




(38) 

[p>+p,(^r^)+p,(r^y+...ydr', 

r'=ro+v. 

Likewise for that part of the disturbed surface which lies below the 
undisturbed surface, equation (10) gives 

F=p / ^inedd I ^^ [Po+Pi(^)+P2(0+ . . .]r^r; (39) 
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For brevity, let the integrals with respect to r in (38) be denoted as 
follows: 

""'- 3{ro+v) ' 






(40) 



•'o = r> > 

Ji'=(ro+A) (ft-c), 

J,'=(r„+t)Moge^, 



(41) 



The ambiguous form which Jl assumes when t=2 receives its proper 
inter])retatiou in the third of (41). 
Similarly, let the integrals with respect to r ia (39) be denoted thus: 

rn (ro+/t)^-ro^ 

"0 = n 1 

Ji"=(»-o+«')A, 

J-/'=(r„+f)Moge!»±^, (42) 

^0 



j-/'=(!g|): rro-('-^)-(ro+/o-(^-^n 



Substituting these equivalents of (40), (41), and (42) in (38) and (39) 
the latter become, respectively. 



>i3 ^27r 



Y=p\ m^eae \ [(Jo+Jo')^o+(«^i+«/i')-Pi+ • • -W^ (43) 
Jo Jo . 

F=p r sin ^^5 r "[ Jo"-Po+ Ji"-Pi + • • • ]^^ . (44) 

(116) 



WOOI.WAKD.] IIAKMONIC DEVELOPMENT OP POTENTIAL. 33 

!Now, it is known from the theory of spherical harmonics^ that the 
general value of P, is 

P,=/ (cos d)fi (cos 6') + terms multiplied by cosines of multiples of (il— A'). 
Since 

/' 'cos i{\^V)d\=% 

we have to deal only with the first term of P<. The function of cos d 
or cos 6' involved in this first term is defined as follows:* 

/i=COS ^, 

The following important relation exists between any three consecu- 
tive values of /, (/i), viz : ^ . 

V V 

The remaining integrals in (43) and (44) are therefore of the form 



Let 



/ (cos 0) Sin eae t ax^^nfi (/) / / {}i)d^ . 

Jo JcoB0 



»^ 008/3 

The known value of this integral is* 

^*(^)=^^[/.i(cos /5)-/,+i(cos/?)] . (47) 

The values of Fin (43) and (44) may now be written thus, replacing 
& by aj (^ or a being the angular distance of the attracted point from 
the center of the attracting mass: 






(/?)] m 



t=0 

for points of disturbed surface above undisturbed. 

%=» 

V=2p7rS^[J/^McoBa)F,(/3)] (49) 

*=o 
for points of disturbed surface below undisturbed. 

^ See Heine, Handbuch der Kugelfanctionen, Theorie und AnwenduDgeD^ Erster 
Bi^nd. Second edition, G. Heimer, Berlin, 1878. 

» Heine, p. 19. » Heine, p. 91. * Heine, p. 93. 

BuU. 48 3 (117) 
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17. The values of Y in the last two equations are exact, but their ap- 
plicability to the problems we have to consider is limited by the slow 
convergence of the series in the second members. In these problems 
^o» ^o+t?, and ro+A are nearly equal, and hence the convergence of the 
series depends almost wholly on the convergence of the functions of a 
and/?. 

If we neglect terms of the order A/ro, t?/ro, and upwards, equations 
(48) and (49) become identical and equivalent to (20) and (21). To show 
these facts we expand J^ and Jl of (40) and (41) and obtain 

j;+c7/=n^[i-4 j (i-l)/^-2t>+(^t+^)|' } +O2], (50) 

in which 0% represents terms of the second order. Likewise, (42) gives 

J-=a[l-^(i-l)^+t 1+0,] . (51) 

Hence, to terms of the first order we have^ 

and (48) and (49) become 



t = QO 



=2ro^p;r \ 

l=;0 



/,(cosa)i^,(/J). (52) 



Kow, if cos ^sscos B cos a+sin Q sin a cos (A— A'), 



t=QO 



..2/.<-«)^.<«=//''vit^ 



■2 cos ^+1 

1=1; " 

But this integral, as shown by the transformation in section 8, is equiva- 
lent to 

/"/cosp-cos/gy X/cogj>-co8^\*^ 

/ vcosj>— cosa/ ^ j vco8i>— cosay 

according as a is less or greater than /?. 

^This inference from (50) and (51) does not appear to be qnite satisfactory. For 
large values of i, Ji '\-Ji' and Ji" are less than roA ; they are each for i=30 . The in* 
fluence of the too great factor ro^ in the higher terms of (52) is, however, counter* 
acted by the small factors /< (cos cil)Fi(J^) in those terms ; and that tbe order of approx* 
jmation secured in (52) is sufficient is evident from the equivalence of (52) with (20) 
and (21), wl^ose order of approximation has been investigated in Article lY, 
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This establishes the equivalence of (52) with (20) and (21), and fur- 
nishes the following development in polar harmonics of the elliptic inte- 
grals in (22) to (25) : 






i=0 



In order to show the form of the second numbers of equations (52) 
and (53) we give below the first four values of / (cos a) and Ft (J3)j re- 
spectively. The series of values may be easily extended by means of 
the relations in (45) and (47). 

/o(co8 ar)=l, Fo{/3)=zl'-cos /?=2 sin^ ^ , 

/i(cos a)=cos a, J'i(/?)= J(l— cos* /?)= J sin* /?, 

/^(cos a)=i{S cos* a-1), J'2(/?)=J(cosy?-co83y?)=J sin*/? cos A 

/3(cos a)=i(6 coV or— 3 cos a)-, JP3(/?)= J(6 cos* /?— 5 cos* /?— 1). 

By means of these values equation (52) may be written thus : 

+2 sin*^+lcos a sin* /3 

2 2 

V=2r^pn |+i(3 cos* «-l) 8in» /? co8 ^ ( ^^^^ 

+■^(5 cos' a— 3 cos a) (6 cos* /?— 6 cos* /?— 1)| 

"P* • • • 

VI. EFFECT OF REARRANaED FREE WATER. 

18. In case the disturbing mass is as large as the supposed ice mass 
of the glacial epoch, the attraction of the rearranged free water on 
itself may be appreciable. To determine the exact effect of this attrac- 
tion would be a work of great difl&culty even if we had the requisite 
information, namely, an accurate knowledge of the complicated shapes 
of the continents and sea Bottom. But we may determine an effect 
which will exceed the probable actual effect by supposing the whole 
surface of the earth covered with a film of water free to assume the 
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proper form for equilibrium under the given forces. To fit this ideal 
case formula (6) may be modified in the following manner: 

Let the potential V in (6) be replaced by F+ JF, where JViB the 
potential due to the rearrangement of the water. Likewise replace 
the constant Vq byFo+^^o. Then, if Jv denote the corresponding 
change in v, equation (6) gives 

v+Av=— — ( F+^ F- Fo- J Fo). (55) 

Now, v+Jv may be expressed by a series of Laplace's functions 
(which are in this case polar harmonics) thus: 

in which Co is a constant of small numerical value ; and hence, denoting 
the deni^ity of sea water by p„ , we have, as shown by Laplace,^ 

^V=WcoP„7t(Zo+lZi+^Z2 ... g-^^i+ . . . )• 

F has already been expressed, equation (54), in a series of Laplace's 
functions. Denoting these functions in (54) for brevity by Fo, Fi, li, 
etc., 

V=2rohp7r{To+Ti+¥2+ . . . ). 

Substituting these values of F, JF, and v+^v in (55) there results, if 
we make the obviously permissible substitution 



2 p„ 4ro/>„;r' 



the following equation : 

r^Zo(l-^)-lh^(To-Uo) 

\, Pm y ^ Pm 

+ rA?,('l-f^-Vi''f K )=0. 

\ *^Pmy ^ Pm 

+ 

According to the theory of Laplace's functions we must have in this 
equation the sums of the functions of the same order separately equal 

^ M^canique Celeste, Book III, Chap. II. 
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to zero, which amounts to placing each line in the equation equal to 
zero. Thus we find 

1*- 

Pm 






{2i+l)p„ 



whence by summation 




Yi 



«^ Pu, 



^^+1 Pm 




(56) 



This equation expresses the total effect of the disturbing mass in 
altering the sea level, v being the effect which would result if the ocean 
were an infinitely rare fluid, and Av being the increase over v which 
would result under the assumed conditions. Obviously, v and Av may 
be expressed separately. Thus 



At:=hB.i V^ ^ Ml^ |. (67) 

Vn. EVALUATION OF CONSTANTS Fq AND Uq. 

19. We proceed now to determine the constants Fo of equation (6) 
and Uq of equations (56) and (57). 

It has already been stated that these constants are to be determined 
from the condition of equality in volumes contained by the disturbed 
and undisturbed surfaces, a condition whose analytical statement is, to 
terms of the order we neglect. 



2rn';r / v sin a(ia=0. 



V si 
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Sabstituting the value of v from equation (6) in this, there results 

/ V sin ador— Fo / sin adar=0. 
t/o t/o 



whence 

sin adCor. (58) 



.=i^r.i, 



(a) The easiest way to evaluate this integral is to substitute for V its 
value given by equation (54). We get, then, at once 

Fo=4roA/);rsin2|, (59) 

since by the theory of spherical harmonics all terms of the series except 
the first.vanish in the integration. For the same reason, if we apply 
the condition 



/ 



{v+/iv) sin a(ia=0 



to equation (56), it will appear that 

I7o=ro=2sin2|. ^ (60) 

(6) The value of F© may also be found by the following process, which 
is chiefly interesting on account of its complication as compared with 
the process used above. For points within the perimeter of the attract- 
ing mass replace Fin (58) by Fi of (20), and for points outside the perim- 
eter replace F by F2 of (21). Making these substitutions, there results 

Fo=2ro&/)f / JiSinada+ / IiB\\iada\. (61) ^ 

Substituting the value of Ji from (24), 




sins^-sins^sinSyi 



l-8in«^8in2yi 



sia*;'i(f ( 8in*0\ 



Pat 

l-sm»^sin»vi 
(122) 
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Then the last integral becomes 



— 4cos*^ 




--4C0S»>^ / r-A_-_A_4.1loff (l+^IOn^n dyr 



in which 



and 



*i=sin|, 



sin ^ cos ^/'i 



("l-sin^lsin^^/iY 

Sabstitating these limits in the non-logarithmic part of the integral, 
it becomes 

48in»|J p—p i^ pf^-mn^^^m^y^'S^dy, 



+2 sin* ^ cot* ^ 



/ sin*;/! "*'«(l-<i)(l+t,v- 



Integrating by parts all terms of this expression except the first, we get 



9 



(l-8in»| sin^x/ ^ \' 

, _^ 2 ^ ^ cot;/-! 

sin -^ sin yi sin ^ 

4 sin'^ ^< + *^^ ®^^ r®^^ ^ ®^^ Ki ) 

— cot^ -^ arc sin (sin ^ sin yi j 



This gives 

Ji sin ada=2(sin /?-/? cos fi). (62) 
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The second integral in (6J ) becomes by sabstitation of the value of 
li from (25) 



r 



I2 sin ada 



p2 

=4 sin^ ^ \ i 

2 J. 



§ \ cos* y2dy2 




<-'0 




(^sin^^^sin*^ sin* yi^Yl^sW iHin^ y^y 



sin ^ cos y2 sin ^ sin* y2 cos ;/2 

ri-sin*|sin*;/2Y Tl-sin*! sinV-jYV 



y2+^ sin 2;/2— 2 arc sin /^sin ~ sin y2 j 
Tl— sin*^ sin*;/2 ) sin y2 



2 



=4sin*^/ sin^ 

2 \ ^ 



arc sin (sin B- sin ;/2 ) 



sin*^ 



2 

This gives 



J2 sin ada=2^7t sin*^— sin /?+/? cos /3\ . ^ (63) 



The sum of (62) and (63) is 



2;r sin* | , 



which, substituted in (61), gives for Vq the same value as (59). 

VIII. EQUATIONS OF DISTURBED SURFACE. 

20. By reference now to equations (3), (6), (20) to (23), and (59) we 
find for the equation of the disturbed surface when the effect of the re- 
arranged water is neglected 

.=3A^§-sia^D. (64) 
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The corresponding expression in polar harmonics is [see equations 
(56), (57), and (60)] 

t=GO 



i=l 



Under the assamption that the water covers the whole sphere and is 
free to adjust itself as stated in Article YI, the equation to the dis- 
turbed surface is 



2 p« 






(66) 



=3AA(£.sin4)+3Ar > ' ' 

^1 L P« J 



21. The position of any point of the disturbed surface is thus defined 
by the co-ordinates v and a^ v being the elevation or depression of the 
point relative to the undisturbed spherical surface and a the angular 
distance of the point from the axis of the disturbing mass. I in (64) 
and (66) is to be computed from (22) or (23) or their equiv;alents (24) 
and (25), accordfbg as the point is within or without the perimeter of 
the disturbing mass. The functions/, (cos a) and ^<(/?) are given by 
(45) and (47), respectively. 

22. The general character of the disturbed surface when the effect of 
the rearranged water is neglected is evident from (64). It is symmetri- 
cal with respect to the axis of the attracting mass. It lies without or 
within the spherical surface of reference according as J is greater or 

less than n sin' ^. The values of I for the point of the disturbed sur- 

face at the center of the attracting mass, for points along the border of 
the mass, and for the point 180^ from the center of the mass, are given 
by equations (26), (28), and (30), respectively. If we denote the corre- 
sponding values of v by the suffixes 1, 2, 3, we get 

.x=3A|-(sin|-8in»§). 
«,=3A A('2 sin» |-siu» |Y 



a=7C. 
(125) 
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tBtJLL.4d. 



The meaning of these equations may be most readily anderstood by- 
reference to Fig. 3. Thus, if the circle FCH represent (in cross- 




FiG. 3. 



section) theMndistarbed sea level surface of the earth, and a stratum of 
matter, as an ice cap, ABDF^ be added thereto, the new sea-level 
surface will assume the form indicated by the dotted line. The values 
of Vi and V2 , as shown in the diagram, are positive, while the value of 
V3 is negative. If, on the other hand, we suppose the space A'B'iyF 
to be occupied by matter of less density than the average density of 
the earth's crust, as is the case in a lake basin, the distui'bed surface 
will fall within the undisturbed surface from F to some line PQ, and 
outside the undisturbed surface from PQ to Jff, i. e., t?i and v^ will be 
negative and Vz positive, or, what amounts to the same thing, p in (67) 
will be essentially negative. 

23. It is of interest to inquire what angular extent of mass will pro- 
duce numerical maxima of i^i , 'V29 and 173, supposing the thickness h and 
the densities p and p^ constant. By means of the usual criteria it is 
readily found that 



i?j=a maximum for /J=60o, 

2 

t72=a maximum for sin y5= -, or /?=39o 32', 

i?3=a maximum for )5=120o. 



(68) 



24. A glance at equation (66) sufSces to show that the effect of the 
free water, if it covers the whole earth, is simply to produce an exag- 
geration of the type of surface defined by (64) and (66). The series in 
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the third member of {^(S) expressiug this exaggeration is rapidly con- 
verging on accoant of the diminishing factor 

3 

2 



(2i+l)^ -3 

Pvf 



which is, since ^ is about — , 

P^ ^ 



22i+5' 



The essential features of the disturbed surface are, therefore, in any 
case, defined by (64) or its equivalent (G5) ^ and in most cases the effect 
of the rearranged water may be neglected as unimportant, or as of no 
greater magnitude than the uncertainties inherent in the data for actual 
problems. 

IX. EVALUATION OF THE DEFINITE INTEGEALS Ji AND I^. 

25. The equations (67) define the position of the disturbed surface in 
some of its most characteristic points. To define its position at any 
other point we must evaluate the elliptic integral Ix or I^^ which per- 
tains to such point. These integrals have already been expressed 
[equation (63)] in a series of polar harmonics, which, if more convergent, 
would suf&ce for computing Ix or I^. It is easy, however, to derive 
more convergent and convenient series than that of (53), and this is 
the object of the present Article. 

First take Ji of (24). For brevity put 

. a 
sm-^ 

sing 



and 



0= sm ^ . 



Then by Maclaurin's series, or by the binomial theorem, we readily find 



(1— -dsin*;/!— Bsin*;/'!— Csin*;/'!— . . •)dyi, 



(69) 



in which 

t '« • 

(127) 
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The even powers of sin yi may each be expanded in a series of the 
form 

c+d cos 2/1+6 cos 4:yi+ , . . 

in which Cy dy e, etc., are constants. Bat since 



/ 



2 

COS 2nydy:=0j 



we shall need in these expansions only the values of c. The valae of o 
in the expansion of (sin yY" is 

_ 2n(2n-l)(2n-2) . . . . jn+l) /1\2» 
^^ 1 . 2 . 3 . . . . n ' \2j 

Applying this formula, and making the integration in (69), there re- 
sults 

Ii=6;r(l-JA-tB-j^eO- . • . ). 

Hence if we put^ 

^i=i(i-n 

^5=^^ai-&^)(7+12&^+186^+286«+636«), 
^g=^j3i^ (1-62) (21+3562+506*+706«+10568+2316W). 



• • • 



Il=•b7r{l—g^ld^—gtW*—g3^tfi— . . . ), 



(70) 



. a 
sin-2 



K?= ^y 6=sin^7 a</3. 



sin§ 

This series converges rapidly, except for values of w near unity. In 
a practical application, to be considered presently, wherein /?=38o (70) 
gives, using terms up to that in k?'^, inclusive, I\ too great by about 5 
per cent, for the case w^\. But this is the most unfavorable case, 

> The general valne of g is 

^4-1 . 1 . 3 . 6 . . . (2n-3) 
+1.1.3.6 . . . (2n-5)l . n6« 

„ _ 2.(a»-l)(2»-2) . ■ . (n+l) J+^-^-^-'' • • • (2»-7)l . 3ft^M 

»"- l..a..8. . . . „. .23» ^+1.1.3.5 . . . (2„-9).1.8.5»V><''-'>y| 

JL . « . o ' 



^-1.3.5 (2n-l)W». 
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and one moreover for which the exact value of Ii is known from equa- 
tion (28). 

26. By a process entirely similar to that followed above, the expan- 
sion of (25) gives, writing for brevity, 

i/=:cosec ^, 

+Vjft«(3+2i^+3»^) 

* +j^^^6«(35+20»'H18i^+20v«+3o»^) 

+T»ttT2&^°(63+35v*+30»^+30K«+35v»+63»'i"')' 

+ 

If in this expression we put 

"'4=2048H"4W?^ +13 1072*^ H" • • • 

245 _L _786 7.2 I 

5= 32768 "T 162144^ "T • • • 

«6= 2 6 2 14 4+ • • • 

we find 

J2=6;r(feitr-*+A:2t(?~'+fe3t<?"*+ . . . ), 

(71) 

K?= ^? 6= sin I , a>/?. 

sin I 

This series converges somewhat more rapidly than (70). For the 
case in which y5=38o and for the extreme value w?=l, using terms to 
that in ti?~" inclusive, (71) gives I2 too small by about 3 per cent. 

•A general expression for the nih term within the brackets, beginning with the 
third term, for which n=2, is the following : 

/+1.3 .5. . .(2n— 1) 
+1.3.5. . . (2ii— 3) 1 . ny2 

^ ^^ i+1 . 3 . 5 . . . (2«— 5) 1.3. ^f.^~~^ ^ y4 
2 ii(2n-l)(2n-2) . . (n+2) V ' ^ 1.2 

l».S?.32...n».2- +1.3.5. . . (2n^7) 1.3.5. ^^ ^^'T^^^^f 



+ 

-f 1 . 3 . 5 . . . (2n— l)K2n 

(129) 
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27. For points near the border of the disturbing mass I2 may be 
expressed by a more rapidly converging series than (71). Thus from 
equation (23) 

l2 



Let 



I \ coa j>— cos a J 
cos /J-cos a=2 sin (^-^ sin (^^r^=2a. 



Ttien 

•^ _ a' a' 




I,= 



fi, ^ ^^_ 3 . "1^ 

I (COS 2?— COS a) (cos p — cos a)^ (cos j) - cos ay I 



Now, if 

r^ dp _ 1 ^_ "^^v. 2 ; 

Jo cosi?— cosa sin« sin f ^~^ \ 



dX '^ I dp 



r 

—sin a I 



da J (cos j>— cos af'^ 



—-_.=— cos a / 7 ^ -+2sin^a / ^ 

dor Jq (cosj)— cosa)* / (cos jp— cos or] 



do^ 

etc.; 
whence 



3^ 






^ dp 1 dX 

(cos jp— cos of sin a da^ 





(iBOSj>— cosar)^ 2sin^aV^^ da)^ 

etc. 

The integrals in the third and higher terms of the above series are 
thus seen to depend on the integral in the second term. Making the 
requisite differentiations we find, to terms of the third order inclusive, 



-/9-/^ ^ .L ^^ cos a ^ a\^-2 sin^ a) ^ \, ^^'^V 2 y 
'"^ ^sETi"*" 2 sin3 ^-t- 4sin«a "*"• ' V ^^'~7«^A 

( 
_/ 5a sin /? , 3a^ cos a sin /? , \ 

V16 8in««+ 8 sin* a +' ' '/ 



a= sin "^ sin ?^, 
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X. SLOPE OF DISTURBED SURFACE. 

28. Having derived the requisite formulas for computiDg the position 
of any point of the disturbed surface, it remains to determine the slope 
of this surface relative to the undisturbed surface. 

Differentiating equation (64) with respect to a^ and dividing the re- 
sult by the radius of the undisturbed surface ro, we get 



dv _ 3hp dl 

r^da^ VfiTtpZ ' da 



(73) 



This expresses the slope or inclination of the disturbed to the undis- 
turbed surface in a meridian plane through the center of the disturbing 
mass ; it also expresses the deflection of the plumb line in the same 
plane. 

In order to apply (73) it is essential to have the general value of 



dl 
da 



Since 



t(?"= 







and hence (70) gives 



3^=-;rcosJA.3^3^5y. (74) 



\ 






Similarly (71) gives 






a~= -^ cos I ^ +f £-3t<»-« ) . (75) 

+ 

29. Equations (74) and (76) will snfBce for the computation of dl/da, 
except for points near to or at the border of the attracting mass. As 
a approaches equality to /3 the above series become less and less con- 
vergent, and finally divergent when a=/3 or w=l. This may be most 
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readily seen by differentiating (22) or (23) with respect to a, and then 
making a—§. Thus we find 







da 2 I cos /^— cos j) 2 



sin ^P"!" 
«'" ^ J. 



Likewise the integrals (24) and (25) become, after differentiating them 
with respect to a and then making a'=y^, 



&a 




dli^ _ I sec^i tanVi(?;/i 

(^sec^|+tan2;/iy' 




da 



^sec^ ^+tan2 ^^ V 



(B) 



IT 




'see* I sec^ Y^iy^ P^ sec^ ;., tan» y^Y^ 



(^sec» l+tanViY / (^8ec» | +tanV»")' 






This shows the equality of (A) and (B) since (B) is plainly infinite, 
its value being 




Xi+rsec'l+tanV^Y 




8ec*4 

It 

30. This failure of equations (74) and (75) for points at the border of 
the attracting mass arises from the fact that the expressions (20) and 
(21), though very approximate for the magnitude of the potential F, 
are not sufficiently general to give an accurate value of dY/da^ or the 
attraction in the direction of the arc a for those points. To determine 
the slope of the disturbed surface at the immediate border of the dis- 
turbing mass a special investigation is requisite. 

Since by equations (3) and (6) the slope is expressed by 

dv __3 1 dV /i^gx 
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we may derive an expression for the attraction d V/roda directly. The 
exact expression for the horizontal attraction towards the axis of the 
mass of any element mass is, asing the same notation as in Article III, 

Ar^dr sin» J cos^ ? dd cos \d\ 

!•' am* I 



r (r— r')*+4rr' sin*-^ j 



and the integral of this is d V/r^da. 
Kowy as heretofore, let 

In addition pat 



so that 



Also let 



whence 



drj^dr^ 
7j=^v for r=ro, 
rf=h'-v for r=ro+ft. 




5=2roSin-, 



dS^ToCOS-^ddy 



\ 



Making these sabstitations and neglecting terms of ttie order 



" * and /^ ^ ^ 



the above expression becomes 

i^dSdrj cos MX 



P 



{S'+rf) 



I 



Integrating with respect to 7, and sabstitutiug the limits given above, 
there results 

r {]i^v)ds , vds 1 ^ ,^ 

Bull. 48.^ 4 (133) 
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If, now, we suppose the attracted point oq the border of the attract- 
ing mass, the limits of S will be 0, and, with sufficient approximation, 
2ro sin p cos il=c cos A, say. Integrating with respect to 5, and substi- 
tuting these limits, we get 

p{h-v) COS UX log. [J 1 + .^j^-^ + -^^-^) 

ui 1 f /TT~?co?T , c COS A\ 
+pv cos XdX loge ^V ^+ — ^^ — + V ) ' 

It remains to integrate these last expressions with respect to A between 
the limits and ~. An application of the formula for integration by 

parts will readily transform them to elliptics, but since their element 
functions decrease very rapidly from the lower to the upper limit, the 
following process^ will suffice. Consider the integral 



Jr^ ^ ji 1 t^ 1 1 . (^ cos^ \ , c cos A^ 
„ cos XdX log, i^yl 1+ -^^- + — ^J 



J 



in which X is such that f — ^—j) ^^y be neglected in comparison with 

unity. In the cases we have to consider ( — ^-^ ) will not exceed -ris 

\ccosXJ 

if cos A=yJ^ or A=89o 25' about. Then, since 



1 r / i . ^ ^ios^ X , c cos A~| 

loe. y 1+ — p- + ^-J 

the above integral becomes 
/ cos XdX loge = loge — I COS XdX + I cos XdX loge cos X 

=( loge — — 1 ) sin A + loge (1+ sin A)+(sin A— 1) loge cos A. 
But since sin A is very nearly unity, the last expression reduces to 

.og.(^)-l. 

Tit 

The error of this integral arising from the use of A instead of -^ as the 
upper limit is less than 

/ ;r ^ \ ^ , / I ^ , c^cos^A , c cos A \ 

which, if cos A=y^ and c cos A/i?=lO, amounts to about FsVTr- 

1 Gfven in a somewhat different form by Helmert in Theorieen der hoheren Geo- 
dasie, Vol. 11, p. 322. 
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For the entire attraction, therefore, of the mass for a p.oint on its 
border we get 

Finally, restoring in this last expression the value of c, viz, c= 
2ro sin fi^ (76) becomes 

f , 8ro sin /? . t? , A— r - 1 ,--, 

I l«&e -kir + A l«g« ^^ 1|- (77) 



dv 3 /tp 



XI. DISTURBED CENTER OF GRAVITY OF EARTH. 

31. Thus far the disturbed surface has been referred to a spherical 
surface concentric with the earth's center of gravity before the dis- 
turbance arose. In determining the effects, of the ice mass in glacial 
times this is the proper surface of reference, since we wish to know the 
distortion of the sea level in those times relative to the sea level in 
preceding and following epochs. If, however, it is desired to consider 
the joint effect in distorting the sea level of existing masses, like the 
continents, on the hypothesis that such masses rest on the surface of a 
centrobaric sphere, a better surface of reference will obviously be the 
disturbed or existing center of gravity of the earth. The use of the 
latter center will require a slight modification of the preceding formulas 
defining the disturbed sea surface. 

To determine the radial displacement of the earth's center of gravity 
due to the addition of such a superficial mass as we have considered, 
it is only necessary to equate the statical moment of that mass to the 
statical moment of the earth's mass, the moment plane being perpen- 
dicular to the axis of the disturbing mass at the undisturbed center of 
gravity of the earth. The moment of an elementary ring of angular 
radius /?', measured froin the axis of the disturbing mass, is to our 
order of approximation 

2u^hp7t sin /J' cos fi'dfi'. 
Hence, if cf denote the displacement sought and M the earth's mass, 

Ma=^n%p7C I 2 sin yS' cos ^'d§' 

=ro%p7r sin^ /3. 

Therefore, by substitution of the value of M given in equation (3), we 
And 

(7=1^-^ sin« /?. (78) 

* Pm 

(136) 
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Now, the elevatioQ of any point of the disturbed surface relative to 
the sphere in the new position will be less than its elevation relative to 
the sphere in the former position by an amount whos* value to the 
proper degree of approximation is 

a cos or, 

a being, as heretofore, the angular distance of the point from the axis 
of the disturbing mass. That is, if i?' denote what v becomes by the 

change in position at the sphere of reference, 

t 

t?'=li— (T COS^. 

Hence, by virtue of (64) and (78) we find for the equation of the dis- 
turbed surface when the sphere of reference is concentric with the dis- 
turbed center of gravity of the earth, 

tj^=3AAr^_sin2 ^("l+cosacos^l-^T (79) 

And the inclination of the disturbed surface to the surface of refer- 
ence is 

^ ^3h-^C^+l sin a sinV\ (80) 

XIL EQUATIONS OF DISTURBED SURFACE WHEN THE DISTUEBING 

MASS IS OF VARIABLE THICKNESS. 

32. Throughout the preceding investigations the thickness of the 
attracting mass has been considered uniform. On this account the 
range of application of the formulas derived is somewhat narrow. It 
may be remarked, however, before proceeding to extend the investiga- 
tion to more complex masses, that inasmuch as the data for actual 
problems will be in general more or less uncertain, or to a large extent 
ideal, formulas of a more comprehensive and hence more complex char- 
acter are not specially desirable. Approximate calculations of a rather 
rough sort in some cases will be as good as the data for those calcula- 
tions. The effects assigned by the foregoing equations will be for the 
most part in excess of the probable actual effects, and in so far as com- 
putation can contribute arguments pertinent to observed facts the max- 
imum effects will be most essential. On the other hand it will be 
desirable in some cases to get an idea of the inferior limiting effects. 
The most important of these cases relates to the extent of submergence 
attributable to the ice cap of the glacial epoch. There would seem to 
be little probability of uniform thickness in such a cap. Apparently 
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some sort of regular decrease ia thickness (with here and there consid- 
erable though comparatively unimportant deviations) from the center 
to the perimeter of the mass is more probable. Such a law of decrease 
is expressed by the equation 



9)(/J)=A=fto 1- 




(81) 



in which h is the thickness along any radial line whose angular distance 
from the axis of the mass is /?, ho is the thickness along the axis, and /3o 
is the angular radius of the perimeter of the mass. In brief, A is a func- 
tion of /3j as stated by the first member of the equation. The exponent 
n must be a positive number and may be for our purposes restricted to 
integer values. In order to determine the effects of masses conform- 
ing to the above, and, in general, any law requiring symmetry of mass 
with respect to a radial axis, we shall devote the present Article to the 
necessary extension of the formulas already derived. 

33. The differential of equation (64) with respect to /3 gives 






dp' 



This expresses the elevation of the disturbed surface due to an annu. 
lus of angular radius /3j of angular width d/?, and height A, the density p 
being uniform. If in this equation we make h a function of /3, or write 
h=(p{/3)y and integrate between the proper limits, the result will be the 
elevation of the disturbed suiface due to a mass whose thickness con- 
forms to the law expressed by q>{/3). Calling, for the sake of distinc- 
tion, v'^ the new value of the elevation of the disturbed surface, and 
the proper limits of /J, /?i, and y52, the result of this integration is 






cp{P)d/3. (83) 



This equation assigns the effect of any homogeneous mass whose 
bounding surface is one of revolution about a radial axis, subject to 
the restriction that the maximum thickness of the mass may be neg- 
lected in comparison with the radius of the earth. It is obvious, 
however, that the integral in (83) may be impracticably complex for 
some forms of q>{/3). To avoid undue complexity and at the same time 
attain results suitable for our special purposes we shall here confine 
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attention to that form of q){^) expressed by equation (81). For this 
function we have, considering the whole mass, 

9>(/J)=0for/?=A, 

n— 1 



sin" ^^ 
2 

Then, observing that J=0 for <3=0, (83) becomes 
^„^3Vr_r. /'Vsin|YVsin|V-A>^sin^4^1. (84> 

■ The definite integral in this expression depends on and will in gen- 
eral be no less complex than J, which is defined by (22) to (26). An 
examination of (24) and (25) shows that for points of the disturbed 
surface within the perimeter of the disturbing mass 



/ J(^sin|y 'dJsin|= / J2(^sin|Y 'eisin 



in /_ 

2 



+ y Ji(^sin|y '(isin|. (85) 



For points of the disturbed surface without the perimeter of the 
mass it is only necessary to replace I in (84) by I2 of (25), or replace 
the limit a in (85) by /Jo. By means of the series (70) and (71), or the 
harmonic series (53), (84) may be evaluated for any point of the dis- 
turbed surface. 

For two points of the disturbed surface, namely, at the center of the 
mass and 180° from that center, (84) yields to direct integration. The 
process of evaluation is as follows : 

The integral in (84) is a function of a, /Jo , and n. Let it be symbolized 
hYf(afton). For the two points noted above let this function be dis- 
tinguished by the suffixes 1 and 2, respectively, so that it becomes 

/i(a, /Jo, n) for a=0, 
and 

/2( «j A> ^) ^^^ OC=7t. 



Then, since from (26), I=7t sin ^ for af=0, 



/i(^,A,n)=-A_;rsin|>. m 
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Likewise, since from (30), I=2:r sin* ^ for a=;r, 



/,(a,/J„»)=2;rsin4"-2^ / ( —7. )dsin»B. (87) 




rsjA -2 COS ^» - /?o cosec ^ ~) for n=l, 
=|A COS* ^« -l\an* ^ for «=2, 



• • • 



When » is more than a few units the integral in (87) may be evalu- 
ated by rapidly converging series. Thus, call the required integral JB, 
and let 



Then 



\ / 



£= 




sin ^ = K sin ^ 



. /3\ » 
8in|- / 










34. It will be particularly essential for our purposes to evaluate (84) 
for points outside the border of the disturbing mass. The integral re- 
quired is, if we write for brevity 

&=sin^ and &o=sin^, 



Now, I2 from the equation preceding (71), may be written thus: 

l2—7r(ji¥+J2¥+jJf+ . . . )i 
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in which ji,J2i etc., are obvious functions of a or i^=cosec a. Therefore 
the above integral becomes 

nC f\b''-'db^n l+i^i^H (89) 

Jo \ -». f 



n 







\T* • • • 

=;r/8, say. 

35. Let the valae of 8 in tbe last equation for points at the border of 
the disturbing mass where a=/?o be denoted by Si . Also denote by 
t7i", 1?/', and i?3" the elevations of the disturbed surface at the center of 
the mass, at its border, and 180^ from its center. Then equation (84), by 
means of the results inl(86) to (89), gives the following equations analo- 
gous to the group (67) : 



^,//=3Ao-^f -^ sin f — ^ sin^ §\ 



a-0; 



V'=3/io^(^.-4^8in«f), (90) 

//=3Ao -^^"2 sin^ ^-2B ^ sin^ &\, 



V3"—6no 



Pn\ 4 n+2 2 f 

a=7r. 



36. To define the slope of the disturbed surface it is in general nec- 
essary to differentiate (84) with respect to a. The result is of a com- 
plex character and subject to discontinuity for points at the border of 
the mass. For practical purposes, however, it will sufiQce to make use 
of Ato"/Aa instead of the differential coefficient, and thus determine 
average slopes over some finite portion of a meridian section of the 
disturbed surface. 

37. In discussing the disturbance of the sea level attributable to tbe 
ice mass of the glacial epoch, it will be of interest to estimate the effect 
of the rearranged free water. For this purpose we may extend tbe 
second term of the third member of equation (^^) so as to make it as- 
sign the effect of the rearranged water when the mass is of variable 
as well as uniform thickness. The process is strictly analogous to that 
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followed in deriving (83) and (84) from (64). Thus, condning attention 
to the form of mass de&ned by (81), and writing, as in section 34^ 



<p{/3)M 1- 




and 



we readily find from (06) 



=Ao(l— 6o""6"), 
JF,(/?)=JF,(6), 



f^'^h.P XI Jo 






1=1 



(2i+l)£--3 

Pw 



(91) 



^ few values of F,(6) derived from (45) and (47) are the following: 

f,(6)=2(6*- h% 

J'3(&)=2(ft«- 66«+10ft«- 56«), 
J'4(6)=2(6»-106*+306«- 356«+ 146«»), . 
^5(6) =2(6' - 156*+ 706«-1406»+ 1266'«-426«). 

From these the corresponding integrals 



/bo 
n— 1 
F,{h)h 



dh 



can be readily derived. Thus, for example, 



+ 



n+2 



V 



nhi 






15 fc4 



n+8 



6o« 



(92) 



+ 



120 . 10 



n+10 
42 



'n+12 



6o^ 
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38. A common property of the formulas (84) to (92), both inclusive, 
is worthy of notice. They all refer to a mass whose thickness conforms 
to the law (81), namely, 

sin«^ 




2 




When »=Q0 this gives (p{/3)=h=hoj or the thickness of the mass is 
uniform. Therefore the formulas (84) to (92) should return to the forms 
applicable to a mass of uniform thickness on making n=x . Such is 
the case. Thus (84) becomes (64), as is readily seen by an application 
of the formula for integration by parts. Likewise, for w=oo equations 
(86), (87), and (89) become (26); (30), and (28), respectively, and the group 
(90) assumes the simpler forms of the group (67). 

B. APPLICATIONS. 

XIII. RELATIVE POSITIONS OP LEVEL OB EQUIPOTENTIAL SURFACES 

IN A LAKE BASIN. 

39. Consider the question stated in section 2, (a), relative to the level 
surfaces in a lake basin. In this case it is required to determine the dif- 
ference in elevation at the center of the basin of two level surfaces 
which intersect along its perimeter. The first two of equations (90) 
give 

i;/'-t?/'=3Ao^f ;Si--A_ sin ^\ (93) 

This represents the difference in elevation of a level or liquid surface 
at the center and at the border of the basin, p must be understood as, 
the excess or defect in density of the liquid relative to the average 
density of the superficial strata of the earth. Thus, if the liquid in 
question be water, 

/)=— (J/}^— 1)=— 1.8, approximately. 

If we differentiate (93), regarding p as variable, the result is 

zr(.."-..")=3ft.f(«:-^siaf). (94) 

This expresses the required separation of the two level surfaces in 
question ; i. e., the separation at the center of the basin of two level sur- 
faces which intersect at the border and which are the free surfaces of 
two liquids whose difference in density is Jp, 

To illustrate more fully the meaning of (93) and (94) let ABDD in 
Fig. 4 represent a cross section through the axis of the basin. ACB is 
a circular arc parallel to the section of the sphere of reference. If the 
basin be filled with water the section of the water surface will lie be. 
low AGB as AG'B. If the water be removed the corresponding sec- 
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tion of the level surface tbrongb A and B will fivll below AG'B as AG"B. 
Hence, iu the diagram 

c 




Fig. 4. 

40. To get some numerical values, let us assume the following data : 

7io=l,000 feet, /?o=arc of lo, A»=5-5, 
p = —1.8, ^P=h 

this latter being approximate difference iu density of air and water. 
Then 

3^0^=982 feet, 

3Ao^=545 feet. 

The assumed value of /3o is equivalent to about 69 miles measured 
along the earth's surface. 

Now, by means of equations (86) to (89), we find the following table 
of results corresponding to several values of n, which defines the shape 
of the basin — see equation (81). The results in the fifth column of the 
table express the difference in elevation of a water surface at the center 
and at the "border of the basin; and those in the sixth column express 
the depression of the level surface at the center of the basin consequent 
upon substituting air for water as the attracting mass. 

Table of values showing relative positions of level surfaces in a lake basin 140 miles in 

diameter and of 1,000 feet maximum (a^al) depth. 



n 


n+1 2 


Si 


n;!''"''-*' 


1 

Vi"-Vi" 


MV2"-'Vl") 


1 


0. 00436 


0. 00161 


0. 00275 


Feet. 
2.70 


Feet. 
L50 


2 


582 


245 


337 


3.31 


1.84 


3 


654 


298 


356 


3.50 


1.94 


4 


608 


333 


365 


3.58 


1.99 


5 


727 


359 


368 


3.61 


2.01 


6 


748 


379 


3G9 


3.62 


2.01 


7 


764 


395 


369 


3.62 


2.01 


8 


776 


407 


369 


3.62 


2.01 


9 


786 


418 


368 


3.61 


2.01 


10 


793 


427 


366 


3.59 


1.99 


00 


0. 00873 


0. 00556 


0. 00317 


3.11 


1.73 
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It will be observed that the differences in the last two columns of the 
above table rise to a maximam between the arguments n=5 and n=9. 
That they should do so is evident from an inspection of equations (89)^ 
(93), and (94). 

The deflection of the plumb line toward the land along the border of 
i uch a lake, supposing it of uniform depth or n=:Qo , would be by 
equation (77) ir' or 16'', according as the basin is filled with water 
or air, and the slope of the water surface at the immediate border 
would be 0.28 feet per mile. 

41. An interesting inference, which might be drawn from the solution, 
of the above problem, is that in triangulating a large lake we should 
expect to see from shore to shore with somewhat less elevations of the 
points of observation than the usual formulas for intervisibility of 
points on the earth's surface would require. This would be a correct 
inference, however, only in case the defect in potential due to the water 
in the lake basin is not offset by an excess in potential due to some 
local or general distribution of matter within the earth's crust. 

XIV. VABIATIONS OF SEA LEVEL ATTRIBUTABLE TO CONTINENTAL 

GLACIERS OR ICE CAPS. 

42. As a second application of the preceding theory we shall investi- 
gate the attractive effects of the ice mass of the glacial epoch, assum- 
ing that the earth's crust did not yield to the pressure of the ice. This 
is the problem of section 2 (b). 

This problem in its physical aspects presents two difficulties, the first 
of which has not been alluded to in the foregoing sections, and the 
second only partially considered. 

The first of these difficulties is to account for the enormous quantity 
of water required to form such an ice mass as is supposed to have cov- 
ered our northern hemisphere during the glacial period. This mass has 
been usually estimated as not less than 5,000 feet thick at its center, 
and to have extended 30° to 90° from that center. It is generally as- 
sumed to have diminished in thickness with some approximation to 
regularity from the center to the perimeter. The superior limit for this 
shape of mass would be a sheet of uniform thickness, and the inferior 
limit a meniscus increasing slowly in thickness from its perimeter 
towards its axis. Evidently, if the water forming such a mass were 
drawn from the ocean, the latter would undergo a considerable diminu- 
tion in elevation, and this diminution might nearly counterbalance the 
attractive effects of the ice in elevating the water along its border. The 
view advanced by Dr. Croll, however, assumes that there is an alterna- 
tion of glaciation at the poles, the epoch of minimum ice cap at the 
one corresponding to the epoch of maximum ice cap at the other. This, 
granting the sufficiency of ice in the two caps, would make the quantity 
of free water in the ocean substantially constant. In the absence of 
definite information on this point, it must be admitted that considerable 
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uncertainty may properly be attributed to our computed variations of 
sea level, although the equivalent lowering of the sea, if the water in 
the ice is drawn therefrom, will be determined for each assumed mass. 

43. The second difSculty arises from the fact that such a large mass 
as we seem compelled to assume for the ice cap would produce an ex- 
tensive rearrangement of the sea water; and we ought therefore, in 
computing the potential at any point of the disturbed surface, to take 
account of this rearrangement. We have shown how to do so in an 
ideal case, which presents effects for the elevation or depression and 
slope of the disturbed surface greater than the probable actual effects. 
The actual effects, as we shall indicate, probably lie about midway be- 
tween those assigned by the formulas for the ideal case and those as- 
signed by the formulas which neglect the potential due to the rearranged 
free water. But on account of the difficulty in fixing an exact limit for 
the actual effects, our computed results will be subject to a sniall range 
of uncertainty, which may be regarded, however, as no greater than 
the inherent uncertainty in the more important data of the problem. 

44. Let us now take for the mean density of the earth p^, for the 
density of ice p, for the thickness of the ice along its axis ho^ and for 
the angular extent of the mass y^oy the following values: 

p^=5.5, ^=10,000 feet, 

p =1,^ /?o=arc of 380. 

This value of /3o corresponds to about 2,600 miles measured along the 
earth's surface ; it is also very nearly that angular extent of mass, of 
uniform thickness, which produces the maximum upheaval of water 
along its border. (See section 23.)^ 

45. As we shall compute the effects of masses corresponding to sev- 
eral values of the index w, equation (81), it will be of interest to define 
with some precision the shape of the exterior surface of each mass. To 
do this it will suffice to give the slope of the surface of any mass at 

^ This value for the density of ice is about 8 per cent, too great ; but by using it 
in the formulas which do not take account of the potential due to the rearranged 
free water, we shall get results differing only slightly from the probable results. *" 

*If the water in the ice cap were drawn wholly from the free sea water the angular 
radius of a mass producing maximum upheaval along its border would be much less 
than the value given by the second of equation {Q8)j for in this case the rise in sea 
level due to the attraction of the cap would be offset partly by the fall due to the 
withdrawal of the water. The difference between the rise and fall ju«t mentioned is 
Vi" of (90), minus r of (96), or 



3*«^(^'-«-T2«''>'t)-ir5:a*«" 



- 2|3o 



The value of (3o which will render this difference a maximum is easily found, but 
the result is of little interest since the corresponding mass would produce effects 
much smaller than the possible effects of continental glaciers. 
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several points in a meridian plane. The differential of equation (81) 
gives 



dh_ 

ndfi 




(95) 



This expresses the slope or inclination in a meridian plane of the 
bounding surface of the attracting mass to the spherical surface of ref- 
erence. Using the above values of ^o and /?©, and for the radius of the 
earth (see section 4) log ro=7.32020, the following table of values has 
been computed. The slopes are expressed in feet per mile. Several of 
the curves whose slopes are given are delineated (with greatly exagge- 
I'ated radial scale) in Fig. 5. 

Table showing meridian slopes of exterior hounding surfaces of assumed ice masses, ISee 

equations (81) and (95).] 



n 






Slopes in feet per mile corresponding to p 


= 




Oo 


50 


10° 


150 


200 


250 


30° 


350 


88° 


1 


Feet. 
3.88 


Feet. 
3.88 


Feet. 
3.86 


Feet. 
3.85 


Feet. 
3.82 


Feet. 
3.78 


Feet. 
3.75 


Feet 
3.70 


Feet 
3.67 


2 


0.00 


1.04 


2.07 


3.08 


4.07 


5.04 


5.96 


6.84 


7.34 


3 


.00 


.21 


.83 


1.85 


3.26 


5.02 


7.10 


9.47 


ILOO 


4 


.00 


.03 


.29 


.99 


2. .32 


4.45 


7.53 


11.66 


14.67 


5 


.00 


.01 


.10 


.49 


1.55 


3.70 


7.48 


13.46 


18.34 


6 


.00 


.00 


.03 


.24 


.99 


2.95 


7.14 


14.92 


22.01 


7 


.00 


.00 


.01 


.11 


.62 


2.29 


6.62 


16. oa 


25.68 


8 


.00 


.00 


.00 


.05 


.47 


1.73 


6.02 


16.97 


29.34 


9 


.00 


.00 


.00 


.02 


.23 


1.30 


5.38 


17.64 


33.01 


10 


.00 


.00 


.00 


.01 


.13 


.96 


4.75 


18.10 


36.68 


00 


0.00 


0.00 


0.00 


0.00 


0.00 


0.00 


0.00 


0.00 


00 



The numbers in the above table show that for w=l the slope of the 
bounding surface, as defined by equation (81), is steepest at the axis of 
the mass and decreases slowly from the axis towards the border. For 
values of n greater than unity the bounding surfaces slope up with 
decreasing rapidity from the border to the axis of the mass, the amount 
of slope diminishing to zero at the axis in each case. The features 
here enumerated will hold for any extent of mass, i. e., for any value of 
j3o. It will be observed also that the slope for any value of n is di- 
rectly proportional to the axial thickness ho» Hence the slopes corre- 
sponding to any other thickness than that assumed (10,000 feet) may be 
readily computed from the table. 

46. As to the actual form of the bounding surface of the ice mass 
of the glacial epoch we have no precise information. It is generally 
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assamed, however, that the mass was thickest along its axis, and that 
the thickness decreased with some approach to regularity between the 
axis and the border.^ The slope of ascent along the border has been 
estimated as 10 to 35 feet per mile.^ Nordenskjold observed a rise of 
7,000 feet in 280 miles from the border of the ice plains of Greenland.^ 




Fig. 5. 

Scale for section of 8pliere=y^7rffoTnTu » 

Scale for distarbiDg mass =^oW=??i)¥oW» 

h=^ho \ 1— ®iHl^^ \ » ;fo=10;000 feet. 
( 8in« iPo ) 

For curve 1, «=1 ; for curve 2, n=2; for curve 3, 71=10 ; for curve 4, n=:x . 

This corresponds to an average slope of 25 feet per mile. It seems 
most probable that the slope of such a mass would be steepest at its 
border and diminish gradually towards its center. Whatever may have 
been the actual slopes it is thought that the preceding table affords a 
suf&ciently comprehensive variety. Our equations, it is true, will as- 
sign the effects of an indefinite variety of other forms, but in the 
absence of more complete actual data the simple forms whose slopes 
have been computed are considered adequate. 

- — III . ■■■■ I ..--_-- — _-■ 

^The directly opposite view is maintained by W J McGee. See his paper on 
Maximum Synchronous Glaciation, Proceedings of the American Association for the 
Advancement of Science, Vol. xxix, 1880. 

* Croll, Climate and Cosmology, p. 244. 

'IWd.,p.245 
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47. In order to form an idea of the amoant of water congealed in tbe 
masses whose shapes have jost been defined we will compute the thick, 
nesses of spherical shells of radios ro (radios of earth's sorface), having 
eqaal volomes with the ice masses, respectively. As the sea covers 
about three-fourths of the earth's surface the products of these thick- 
nesses by ^ will represent approximately the necessary lowering of the 
sea level if the water in the ice cap is drawn from the sea. Since the 
thickness of our assumed ice mass along any radial line at an angular 
distance fi from its axis is by equation (81) 



h:=z 




sm 




the volume of this mass will be expressed by the integral 



2ro^ho7r 




sin /3d/3=4: -?_ n^hn: sin» A. 



The thickness r of a spherical shell of radius n and equal volume 
with the above is given to a sufficient degree of approximation by the 
relation 



whence 



Wr7r=4: _?- ro'hoTt sin^ ^ , 



r=; 



n 



n+2 2 



(96) 



From this equation with our working values, ^0=10,000 feet and /3o= 
380, we find the following values of r and |r, corresponding to tbe 
values of n in the preceding table : 

TdbU showing thicknesses r of spherical shells of equal volume with assumed ice masses, and 

equivalent lowering of sea level f r. 



n 


T 


■|T 




Feet. 


Feet. 


1 


353 


471 


2 


530 


707 


3 


«36 


848 


4 


707 


943 


5 


757 


1009 


6 


795 


1060 


7 


824 


1099 


8 


848 


1131 


9 


867 


1156 


10 


883 


11T7 


oo 


1060 


1413 
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48. Although, as seen from the last table, the amount of water nee* 
essary to form such an ice capas we are considering would, if drawn 
from the sea, cause a decided lowering thereof, yet the mass of ice 
would be very small compared with the mass of the earth. Thus, for 
example, the mass of a sheet of uniform thickness 10,000 feet and 38<> 
angular radius is only ye^tnr P^rt of the earth's mass. 

49. To determine the position of the disturbed relative to the undis- 
turbed surface it will be suflBcient to compute the elevations of the 
water at the center of the attracting mass, along its border, and 180<^ 
from the center, by means of formulas (90). As the slope of the dis- 
turbed surface near the border of the mass is of most importance, we 
shall compute the elevation of the disturbed surface for a circle of 
points 1° distant from the border, and thereby deduce the average 
slope of the disturbed surface within that distance (69 miles) of the 
border. The separate quantities required in this calculation are given 
in the following table for the same values of n as those used in the two 
preceding tables. It will be remembered that ti=QO corresponds to an 
attracting mass of uniform thickness. The values of S in the last col- 
umn of the table have all been computed from equation (89), except the 
one for n=cc , which has been derived from (72). 

Table of numerical values of functions defining position of disturbed surface^ 



n 


n+1 2 


n+2 2 


Si 


1 


0. 16278 


0.03533 


0. 06078 


2 


.21704 


.05300 


.09249 


3 


.24418 


.06360 


.11218 


4 


: 26045 


.07066 


.12568 


5 


.27131 


.07571 


. 13555 


6 


.27906 


.07950 


. 14310 


7 


.28487 


.08244 


. 14908 


8 


.28939 


.08480 


.15393 





.29301 


. 08672 


. 15795 


10 


.29597 


.08833 


. 16135 


00 


0.32557 


0. 10600 


0.21111 



2(8in»^-B) 


8 


0. 01796 




0.05885 


. 02699 




.08946 


. 03243 




. 10842 


.03606 




M2139 


. 03867 




.13085 


. 04063 




.13809 


. 04216 




.14380 


. 04338 




.14843 


. 04438 




.15227 


.04522 




.15550 


0.05448 




0. 19846 



^ For brevity make the following substitutions: 



^ n-fl 2 n+2 2 ^ 



iV2=/8i- 



n 



n+2 



sin^^ 



A 



TT ' 



S'3=2('8m^^-B\ 



n 



n+2 



sin* & 



2"' 



N^=8i-S. 



Ball. 48- 



(149) 



66 



FORM AND POSITION OF THE SEA LEVEL. 



[BULL. 48. 



TheD from the preceding table we get the following table of values: 



n 


JV^ 


.V, 


^Va 


X* 


1 


4-0. 12745 


+0. 02545 


-0. 01737 


1 

4-0. 00193 


2 


+ .16404 


4- .03949 


- .0:^601 


4- .00303 


3 


+ .18058 


+ .04858 


- .03117 


4- .00376 


4 


+ . 18079 


+ .05502 


- .03460 


4- .00429 


5 


+ .10560 


+ .05984 


- .03704 


4- .00470 


6 


+ . 19956 


4- .16360 


- .03887 


+ .ooroi 


7 


+ .20243 


4- .0e664 


- .04028 


+ .00528 


8 


+ . 204:.9 


4- .06913 


- .01142 


4- .00550 


9 


+ .20629 


4- .07123 


- .04234 


4- .C0568 


10 


+ .20764 


4- .07302 


- .04311 


4- .00585 


00 


+0.21957 


4-0.10511 


-0. 05152 


4-0. 01265 



Now, the factor by which we must multiply Ni, N2J and N3 to get the 
elevation of the disturbed surface above the undisturbed at the center 
of the mass, along its border, and 180^ from its center is 

3Ao^ =5454.5 feet: 

and one sixty-ninth part of this factor multiplied by .5^4 will give the aver- 
age slope per mile of the disturbed surface within 1^ of the border of 
the ice. Hence we get in the table below the results corresponding 
to the several values of n. The plus sign indicates elevation and the 
minus sign depression of the disturbed relative to the undisturbed 
surface. 

Table showing effects in distorting the sea level of ice caps of the same angular radius j 38°, 
and same axial thickness, 10,000 feet, hut of varying external slopes j defined by equatidn (81). 



n 


Position of disturbed relative to undis- 
turbed surface. 


Average 

slope per 
mile of dis- 
turbed sur- 
face within 

10 of the 
border of the 

ice mass. 


At centerof 
ice mas. 


A Ion 2 border 
of ice mass. 


180= from cen- 
ter of ice 
mass. 




Feet. 


Feet. 


Feet. 


Feet. 


1 


4- 605 


+139 


— 95 


0.15 


2 


4- 895 


+215 


—142 


.24 


3 


+ 985 


+265 


—170 


.30 


4 


4-1035 


+300 


—189 


.34 


5 


+ 1067 


+326 


—202 


.37 


6 


4-1088 


+347 


—212 


.40 


7 


4-1104 


+363 


—220 


.42 


8 


4-1116 


+377 


—226 


.43 


9 


4-1125 


+389 


—231 


.45 


10 


+1133 


+398 


-235 


.46 


00 


+1198 


+573 


—281 


1.00 
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50. As already explained, the numbers in the second, tbird, and fourth 
columns of the last table assign the position, and those m the fifth col- 
umn the meridian slope of the disturbed surface relative to the undis- 
turbed surface, assuming that the ratio of the density of ice to the earth's 
mean density is -^^f , and neglecting the effect of the rearranged free 
water. If ^ were the correct ratio of the densities, and if the sea cov- 
ered the whole surface of the earth, and were free to arrange itself in con- 
form ity with tUie attractive forces, formulas (66) and (91) show that the 
numbers in the second, third, fourth, and fifth columns of the table 
should be increased by about 15, 18, 26, afid 18 per cent., respectively, of 
their stated amounts. But the ratio -^^ is too great by about 8 per cent., 
so that for this reason the above percentages must be reduced to 7, 10, 
18, and 10, respectively. Again, only three-fourths, at most, of the 
earth's surface is covered with water, so that the reduced percentages 
must be diminished to three fourths their stated amounts. Moreover, 
the sea could not penetrate the ice mass in such a manner as to produce 
the assumed increase in the potential within its border, and hence the 
reduced percentages must be still further diminished. It is estimated 
that the actual effect can not be greater than two-thirds that which 
would follow if the water were unrestricted. Accordingly, § . f = J the 
above reduced percentages, or 3.5, 5, 9, and 5 per cent., respectively, 
would appear to be liberal allowances for the effect of the free water 
in exaggerating the deviation of the disturbed from the undisturbed 
surface at the points designated, over and above the tabular devia- 
tions. In view of the smallness of these possible increments to our 
computed quantities we need give the question of the effect of the re- 
arranged free water no further consideration. 

The results in the second column of the table are the heights to 
which the water would rise at the center of the ice mass if brought 
within it in any manner, as by a canal, and left free to assume equi- 
librium. If the amount of free water were sufficient it would rise or 
fall to the extent indicated in the third and fourth columns, respectively, 
at points, along the border and at the antipodes of the center of the 
mass. 

The slopes given in the fifth column will apply to isolated bodies of 
free water adjacent to the ice mass, and also to the sea surface in the 
same vicinity, whether there be sufficient water to rise to the height 
indicated in the second column or not.^ 



^Recent observations of the beaches of bodies of water contiguous with the ice 
fields of the glacial epoch indicate uniformly that the water surfaces sloped upwards 
toward the ice. In a careful exploration of the beaches of Lake Agassiz, an extinct 
lake which lay in the valley of the Red River of the North during the glacial epoch, 
Professor Upham has found slopes varying from zero to 1.3 feet per mile. (See Bulle- 
tin No. 39 of the U. S. Geological Survey, on The Upper Beaches and Deltas of the 
Glacial Lake Agassiz, by Warren Upham.) Near the south shore of Lake Ontario, in 
New York, Mr. G. K. Gilbert has observed slopes as great as 5 feet per mile. (Sea 
Science, Vol. I, p. 222. ) 
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The maxiinam possible slope of tbe disturbed surface would occur at 
the immediate border of a mass of uniform thickuess. This slope is by 
equation (77) for ^=10,000 feet and /J=38o, 1.80 feet per mile, and cor- 
responds to a plumb line deflection of 72''. Since it is not probable 
that the ice cap presented at its border anything like a vertical wall 
10,000 feet high, we infer that a mass whose maximum or axial thick- 
ness is 10,000 feet would be quite inadequate to produce a slope of 1.8 
feet per mile. 

Again, since all the results in the table are proportional to the axial 
thickness of the ice, to produce as great an average slope as 2 ieet 
per mile within 1^ (69 miles) of the border of a mass having the more 
probable slope defined by the index w=6 to w=10, would require an 
axial thickness in round numbers of 50,000 feet, or ^} miles. The 
slope at the immediate border in the extremely improbable ca^e of a 
uniform thickness of 60,000 feet would be 5x1.8 feet =9 feet per mile. 

51. To produce slopes as great as 4 or 5 feet per mile for any distance 
from its immediate border an ice cap must have great thickness, which 
implies for any large areal extent a heavy draft on the visible supply 
of water. The minimum thicknesses of ice masses of varying radial 
extent which would produce an average slope of 5 feet per mile within 
lo (69 miles) of their borders are shown in the following table. These 
values are computed on the improbable supposition that the masses 
are of uniform thickness. The volume of each mass is indicated by the 
equivalent lowering of the sea level, or |r deduced from equation (96), 
the value of n being infinite. 

Table showing for ice masses of varying radial extent the minimum thicknesses requisite to 
produce an average slope of 5 feet per mile within 1° of the borders. 



Angnlar radins 
of mass. 


Minimnm thick- 
ness. 


Equivalent 
lowering of sea 
level, indicat- 

inn volume 
of masa. 


o 


Feet. 


Feet. 


10 


69, 400 


701 


20 


52,500 


2,3U8 


30 


52, 600 


. 4, 699 


38 


50, COO 


7,065 



For masses having moderate surface slopes near their perimeters the 

axial thicknesses must be about twice as great as the minimum values 

given in the table to produce the same average slope of 6 feet per mile 

within lo of the borders. Thus, for the angular extent 38^, and 

for w=6, say, an average slope of 5 feet per mile would require an 

axial thickness of about 125,000 feet, or 24 miles. This corresponds ta 

a lowering o/ the sea of about 2.4 miles, and if the quantity of free 

water were suflBcient, to an elevation of the sea along the border of the 

mass of about 4,000 feet. 
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52. Confiniug further investigation to our working axial thickness, 
10,000 feet, angular radius 38o, and ratio of densities t^, we may in- 
quire as to the extent of the variation in sea level at any point of the 
earth's surface on the supposition of an alternation of glaciation at the 
poles. For this purpose it is simply necessary to compute the eleva- 
tion of the disturbed surface at the angular distance a of the point in 
question from the pole or axis of the ice cap and for the point 180O— a 
and take the difference between the results. Since the maximum effects 
are of most interest, we shall compute the variations in sea level on the 
assumption that the ice cap is of uniform thickness. The results, as 
may be inferred from the second table on p. 6Gj will not differ materi- 
ally, except for points near the border of the ice, from the results which 
would be derived on the assumption of a sloping mass corresponding 
to the index n=6 to n=10. 

To compute the required integrals we may use formulas (70) and (71). 
For /3=38o we find 

log ^1=9.3493, log A;i=:9.7049, 

log ^2=8.7422, log A'2=8.8079, 

log ^3=8.3452, log A'3=8.3849, 

log<74=8.0708, log ^'4=8.1186, 

log ^5=7.8644, log A*5=7.8904, 

log ^6=7.6987. log A'6=7.7033. 

From the formula 

sin I 

we compute log w for a=lO^, 20^, 30^, etc., to 180^ and then the values 
J, or, more conveniently, the values of I/tt sin c required in (64), readily 

follow. The results are grven in tabular form below for circles at inter- 
vals of 10^ from either pole. 
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Table showing dislitrhance of sea level attributable to an ice cap of 38^ angular radius and 
lOjOOO feet uniform thickness, and variation in sea level attributable to same nmss on the 
hypothesis of an alternation of glaciation at the twopotes. 



Angular distance 

from either pole, 

or a. 



Elevation of dis- T*"? 1*^ 'u ^* ^^^-^ ! 

lurbed above nn- i^^'",^P***'*'^^°"'**- 

disturbed surface. "^^'^ f pooh of max- . 
-_ ^ ' imum glaciation, or 

la — cigo — a. 



o 


10 
20 



Feet. 

+ lll,»8 

+ 1169 
-fl076 



Feet. 
Ii79 

1449 

1352 



30 


+ 893 


iiaj 


40 


+ 441 


702 


50 


+ 185 


4:u 


60 


+ 46 


279 


70 


- 43 


170 


80 


- Il'6 


80 


90 


- i:.2 


iO 


lOJ 


- 186 


80 


no 


— 213 


170 


120 


- 233 


279 


130 


- 249 


434 


140 


- 2GI 


702 


150 


— 270 


1163 


160 


- 276 


1352 


170 


- 280 


1419 


180 


- 281 


1479 



A graphical represeutatioa of the resalts in the last two coluinus of 
the above table is given in Fig. 6, the section of the disturbed and un- 
disturbed surfaces being so developed that the great circle of the latter 
appears as a straight line, AS. The distances between the two curves 
representing the disturbed surfaces measured at right angles to the 
axis AB indicate the variation in sea level between the epochs of max- 
imum and minimum glaciation at either pole. 

The variation in slope of the sea surface at an3^ point during the in- 
terval between the extremes of glaciation will equal the sum of the 
slopes at that point for the two epochs. The maximum variation is, 
however, only slightly greater (about 1 per cent.) than the maximum 
slope already computed, and requires, therefore, no further considera- 
tion. That portion of the disturbed surface having the greatest slope 
is made clearly apparent by the curves in Fig. 6. 




JhsturbeeL 



Surface. 



Fig. 6. 
Horizontal scale, Tsmhwaji^^ vertical scale, TTrfSSSSoir* 
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XV. HISTORICAL NOTE. 

53. The effect of the glacial accumulation ia disturbing the sea level 
Las been the subject gf considerable discussion. The discussion was 
started apparently by Dr. GroU about 1866; and from the widely differ- 
ing results quoted by him in his Climate and Time, Chapter XXIV, one 
might infer that the mathematicians who have attacked the problem 
are completely at loggerheads. There has been, indeed, some diversity 
of opinion as to the proper method of treating the question, and the 
work of one writer at least is quite erroneous; but part of the appar- 
ent discrepancy in the results quoted by Dr. Croll is due also to a rad- 
ical difference in tlie data used as a basis for calculation. 

The mathematicians whose writings on this subject appear to be 
worthy of especial consideration are Archdeacon Pratt, Mr. D. D. 
Heath, and Sir William Thomson. Their investigations may be found 
in the Philosophical Magazine for 1866 ; Vols. XXXI and XXXII. 
^hose of Archdeacon Pratt are reproduced substantially in his Figure 
of the Earth, fourth edition, pages 236 to 238, while those of Sir Will- 
iam Thomson are given also in Croll's Climate and Time, Chapter XXIII. 

54. Pratt — In the following notice of Pratt's investigations we shall 
refer to the fourth edition of The Figure of the Earth (published by 
McMillan & Co., London and 'New York, 1871). as it is the more recent 
and formal treatise. Although in many respects a valuable text book, 
it is marred by some serious errors, one of which we proceed to point out. 

Pratt's investigation is based on the foil owing erroneous proposition 
(see page 212, Figure of the Earth), namely : 

To prove that the effect of a mass at the earth's surface, whether above or below, 
is to make the sea level rise at any place through a space, , where Fis the poten- 
tial of the mass for a point on the disturbed sea level, which is in the same vertical 
line with the place. 

By a process of reasoning to which there appears to be no objection 

he arrives at this equation, namely, 

Y 
r+coust= — , (a) 

9 

in which r is the radius-vector of any point on the disturbed sea sur- 
face, y the potential at that point of the disturbing mass, and g is the 
well known velocity increment due to the earth's attraction. The radius 
of the undisturbed surface, supposed spherical, being denoted by a, 
Pratt says : 

Let r=a where F=0 or the horizontal attraction of *the mass first becomes appre- 
ciable. 

He thus finds in the above equation 

const= — a, 
and hence 

rise in sea level=r— «:^ -- . (h) 

g ^ 
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But this reasoning is strangely faalty. V is not zero for r=a j in 
fact it is zero only for points infinitely removed from the surface under 
consideration. He should have reasoned thus : Along the line of inter- 
section of the disturbed and undisturbed surfaces r=a. Call the par- 
ticular value of F along this line Vq. Then we have from (a) 



y 

const=— a+ —. 

9 



and hence 



9 

This agrees with our formula (6)', and the constant Vq is to be deter- 
mined from the condition that the disturbed and undisturbed surfaces 
contain equal volumes. Formula (a) gives results which are too great 

by the constant amount — , i. e., this formula measures heights above 

y ^ 
a spherical surface — below the undisturbed surface (see section 5). 

"To, it will be observed, is never of a lower order than V and can not 
therefore be neglected in comparison with V. 

Although Pratt makes use of the correct principle for determining 
the constant Fo in his article 199, he ignores this principle altogether 
in his article 200, referred to above, and again in his article 213. 

55. In his calculation Pratt assumes a sheet of ice 7,000 feet thick at 
the pole to extend over a whole hemisphere, decreasing in thickness^ 
however, as it recedes from the pole in the ratio of the square of the 
cosine of the polar distance. He does not consider the effect of the 
rearranged water. His method of determining the potential V is not sat- 
isfactory. It consists (see articles 90, 91, and 92, Figure of the Earth) 
in a species of mechanical quadrature, by which he computes five special 
values of the attraction of a '* hemisphericospheroidal meniscus," whose 
thickness varies according to the law stated above. From these five 
values he derives by the method of indeterminate coefficients a general 
formula for the attraction of the meniscus, and from this formula by 
integration he gets a general formula for the potential of the mass. 
The order of approximation of these formulas is not shown and is not 
evident. From a test we shall ax)ply it is inferred that the approxima- 
tion is so rough as to rend^er the formulas worthless. 
. We may readily derive the proper expression for the elevation of the 
sea under the conditions assumed by Pratt from our general equation 
(83). In this equation, if the thickness at the pole be denoted by ^o? we 
have 

(loO) 
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Therefore (83) becomes 

^Pn^Jo dp ^ ^ 2 p^ (d) 

Y 
The first term in the second member of this equation is — , and the 

second term is — of equation {c) above. Ifo w, for the ratio of the density 

of ice to the mean density of the earth ^ Pratt uses - , and hence the 

Pm ^ 

constant by which his results shouUI be diminished, if they were what 
they purport to be, is 

^^7000 feet=583 feet. 
12 12 

xr 

56. We may test the correctness of Pratt's formula for computing — 

by deriving the elevation of the disturbed surface at the center of his 
ice sheet. For this point we have by equation (26) 

2 ' 
and hence 

___cos^. 
The integral in (d), therefore, is 

5/^03 |e»s.^,=.[Q)*-|(J)'+J(>)'] =4, g.; 



and (d) becomes 



*"=ir'^»-r2'^»- ^^y 



This assigns the height of the disturbed above the undisturbed surface 
at the center of the ice mass. Now, 

J^ =0.18856, 
li> 
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but Pratt's formula (see page 237, Figure of the Earth) gives 0.1189, which 
is only about 63 per cent, of the true value.^ By a fortuitous compen- 
sation of errors, however, Pratt's formula gives a nearly correct result 
for the elevation of the sea at the center of the ice cap, the error in the 
potential for this point being about eqaal to the constant omitted. His 
formula gives for this elevation 0.1189/to=0.1189x7000 feet=832 feet. 
The correct value is (0.18856—^)^0=737 feet, which is about 13 per 
cent, of itself smaller than Pratt's value. 

57. Heath. — The investigations of Heath may be found in numbers 
covin, COIX, and CCXIII of the Philosophical Magazine. In num- 
ber covin he develops the theory of the effect on the sea level 
of an ice cap of uniform thickness and any angular radial extent, 
and applies this theory to a numerical example. In number COIX he 
corrects a blunder by which he was led in his previous paper to the 
conclusion that the ice cap would produce a rise of sea level in both 
hemispheres. His last paper in number CCXIII is chiefly interesting 
as a review* and criticism of the work of Oroll, Archdeacon Pratt, and 
Sir William Thomson. 

The method followed by Heath is that of Laplace's functions (spher- 
ical harmonics). He takes account of the rearranged water on the 
supposition that it covers the whole sphere and is free to permeate the 
ice mass. He considers masses of uniform thickness only, arguing, 
though not very cogently, that this was the probable form of the ice 
cap. His mathematical processes are correct in principle, but his for- 
mulas defining the position of any point of the disturbed surface are 
rather uninviting to the computer on account of the slow convergence 
of the series used. His series is the same as that in the second mem- 
ber of our equation (66), which we have separated into the definite in- 
tegral and rapidly converging series of the third member of (66). 

iThat the first term in the second member of {d)' is — for the point in question is 


easily proved directly. Thus, to terms of the order we neglect, the potential of the 

mass is 






Y^ I ^ 2nro sin /3 . ^p cos^ /? . rpdp 

2ro sin ^ 






15 
Dividing this by <7=}ropmn- we get 



=2 )'ohopiT Iq cos 2 cos2 j3dlS 

^4/'o«/J7r- 
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58. In his numerical example Heath assumes an ice sheet of uniform 
thickness, 10,000 feet, to extend 30^ in every direction from the pole. 
For the ratio of the density of ice or sea water to the mean density of 
the earth he uses i. That is, in the notation .of our equation (66) he 
has 

/t= 10,000 feet, 

/?=30o, 

Heath gives the position of the disturbed surface at the two poles 
and at points 35^ distant from either pole. He uses fifty terms of his 
series. We will verify one of his results, viz, that assigning the eleva- 
tion of the disturbed surface at the center of the ice cap. With the 
above data formula (66) becomes 

^+Jt>=g^(j-;rsin^|)+250oV(/'(««».y'(/^)). 



l-.l 



At the center of the ice cap 



a'=0, 

J=7r sin ^, 



/i (cos 0')=1. 
The numerical values of the first ten terms of 






are, in order. 



+0.02500, 
+ .01203, 
+ .00661, 
+ .00358, 
+ .00172, 



+0.00058, 

- .00008, 

- .00039, 

- .00047, 

- .00041. 



The sum of these is 



+0.04817. 



Therefore the above equation becomes 

r+z/t?=959 feet +120 feet = 1,079 feet. 

For the same result Mr. Heath gives 1,078 feet. 

69. In his review Heath first devotes some space to a criticism of the 
views of OroU and the work of Thomson. CroU had apparently held 
the notion that the change in sea level due to an ice cap is essentially 

(159) 
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equal to the change in position of the center of gravity of the earth 
and cap. This notion is indeed correct for the particular ideal case con- 
sidered by Croll and Thomson, but it is in general quite incorrect. 
Henth points out this fact and concludes that the special case discussed 
by Croll and Thomson affords us an inadequate conception of the actual 
problem. 

With reference to the investigation of Archdeacon Pratt, which hi 
been published in number CCVIII of the magazine, Heath says : 

I mast confess, with some diffidence, that it appears to me radicaUy erroneoos. 

He then proceeds to mention, correctly in the main, but with some 
reservation, the defects of Pratt's method. On one point, however, he 
expresses some doubt as to the adequacy of his own process. He is 
not sure that the order of approximation of his expression for the po- 
tential of the ice cap is sufficient, and thinks this may require closer 
investigation. We have cleared up this point in Article IV. 

60. Thomson. — The contribution of Sir William Thomson to the dis- 
cussion of this subject is appended in the form of a note to a paper by 
Croll, On the Physical Cause of the Submergence and Emergence of 
the Laud During the Glacial Epoch, published in number CCIX of the 
Philosophical Magazine. The note is brief, but it contains a clear state- 
ment of the essential analytical considerations required in the solution 
of the general problem, and of the decided simplification which results 
when the ice sheet has the special form assumed for the purpose of dis- 
cussion by Croll, namely, that of a hemispherical meniscus, whose thick- 
ness (or density) varies everywhere as the sine of the latitude. In 
analo;:jy with the views of Croll, Thomson devises the following ideal 
conditions, which are interesting as presenting the mechanical features 
of the problem in their simplest form. He says: 

As an assumption leading to a simple calculation, let ns suppose the solid earth to 
rise out of the water in a vast number of small flat- topped Islands, each bounded by 
A perpendicular cliff, and let the proportion of the water area to the whole be eqnal 
in all parts. Let all of these islands in one hemisphere be covered with ice, of thick- 
ness according to the law assumed by Mr. Croll, that is, varying in simple proportion 
of the sine of the latitude. Let this ice be removed from the first hemisphere and 
similarly distributed over the islands in the second. 

Thomson gives no -analysis, but continues: 

By working out according to Mr. CrolVs directions, it is easily found that the 
change in sea level which this will produce will consist in a sinking in the first hemi- 
sphere and rising in the second through heights varying according to the same law 
(that is, simple proportionality to sines of latitudes), and amounting at each pole to 

(l-GO)U 

1— (»«7 ' 

where * denotes the thickness of the ice crust at the pole, i the ratio of the density 
of ice, and w that of sea water to the earth's mean density, and go the ratio of the 
area of ocean to the whole surface. 

61. We may readily get Thomson's result from equation (66). Thus, 
4 the change in position of any point of the sea surface will be, if we 

^ (16:)) 



WOODWARD.] INVESTIGATIONS OF THOMSON. 77 

represent the tliickuess of the ice at auy point whose polar distance is 
/? by *o cos /?, 



/ 



^(^A.C08/Se7/J= 



^ //.(cos a) r ^^) cos fid^ 



2**^ /j\ 1_ 3 p, 




Bat from equation (46) 

^W)=y;(cosy8)sin/?; 

and by the theory of Laplace's functions 



/ 



f, (COS /?) COS /? sin /?(i/?=0, 





except when i=l. In this case, since /^(cos y?)=cos /?, the last integral 
becomes 



/ 



IT 

coB^/3 sin /3d/3=^. 

^ 



Therefore, observing that/ (cos a)=cos a, a being the polar distance 
of any point of the sea surface, we. get 



/ 



^ ^Q ho cos /3d/3=: -t^ — cos a. (e) 



'to 



'm 



Now, since the ratio of the area of the ocean to the whole surface of 
the earth is assumed in Thomson's problem to be c«7, we must replace /> 
in {e) by (1— c»)p and p„ by gdp„. Making these substitutions, and put- 
ting a=0, the second member of (e) becomes 

^- (/) 

This is Thomson's result. 

• Knowing the fact expressed by (e), namely, that the transfer of such 

a meniscus as we are considering from one hemisphere to the opposite 

one would change the sea level at any point by an amount proporf iouail 

to the cosine of the polar distance of that point, we may get^the result 

(161) 
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(/) ^y equating the sum of the moments of the transferred ice and 
water to the moment of the whole earth. In other words, this partic- 
ular case requires only the amount of shifting of the earth's center of 
gravity. Calling this amount (Ti , and taking a plane perpendicular to 
the axis of the meniscus at the disturbed center of gravity of the earth 
as moment plane, the equation of moments is 



(l-a,)pj; 



2;rro sin 6 . ndd . //q cos 6 . ro cos 6 + 



Gt)f)„ I 27rro sin 6 , rdH . Ci cos 6 . r cos 0=-o ^^oVw^i* 



This gives 

0-1 = -^, 



P 



M 



which is the same as (/). 
In his numerical example Thomson takes 

/io=G,000 feet, 

2 

^4, and A. 2^ 

Pm 6' p^ 11 

These data give 379 feet as the change in sea level at the pole during 
the interval between the epochs of minimum and maximum glaciation. 
This is the greatest change in sea level that could occur under the as- 
sumed conditions. 

62. It is to be observed that the numerical results of Pratt(correcte(l), 
Heath, and Thomson are not directly comparable with each other nor 
with the results we have computed in Article XIV, since they are all 
based on different data. They represent effects due to causes of the 
same kind, but of widely differing magnitudes. No statement of these 
results would be intelligible without an accompanying statement of the 
data on which they rest. To show clearly how widely different writers 
differ in their data the latter have been collected in a tabular form be- 
low. It will be remembered that Pratt took no account of the effect of 
the rearranged water :^ 

1 One of the most important contributions to the disoassion of the effects of conti- 
nental glaciers on the sea level has appeared since the mauasoript of this paper was 
placed in the printer's hands, viz : Die Geoiddeformationen der Eiszeit, von Erich 
von Drygalski, Dr. Phil. W. Pormetter, Berlin, 1887. 
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Table showing data used by different authors i?i discussing the problem of glacial submer- 
gence. 



Author. 

1 


Ratio of 

density of 

ice to mean 

density of 

earth. 

P^ 
P.^ 


Ratio of 

density of 

sea water to 

mean density 

of earth. 

pa 
Pm 


Angular 
radius of 
ice mass. 


Thiokne«8 of ice at 

angular distance /3 from 

its axis. 


Pratt 


h 
h 
i 






00 
30 
90 

38 


Feet. 
7,000 co8*/3 
10, 060 uniform. 
6, 000 cos /3. 


! Heath 

Thomson 

Woodward 


h 
ft 



XYI. VARIATIONS IN SEA LEVEL ATTRIBUTABLE TO CONTINENTAL 

MASSES. 

63. To illustrate the use of the theory developed in determining the 
disturbance of the sea level attributable to a continental mass we shall 
consider, in addition to some observations on the general features of 
the subject, the special case presented by the largest of the continents, 
namely, Europe and Asia. 

In contemplating this problem it is important to distinguish two ex- 
treme hypotheses relative to the nature of the earth's crust. On the 
one hand, we may suppose that a continent is simply a superficial ag- 
gregation of matter, which, if removed, would leave a sensibly centro- 
baric spheroid } in other words, the presence of a continent does not 
imply that beneath it the earth's crust is any less dense than beneath 
the ocean. On the other hand, we may suppose that the several radial 
element prisms of the earth's crust are in a state approximating to hy- 
drostatic equilibrium, and hence, that the mere existence of a continent 
implies a defect of density in the strata beneath it. According as we 
proceed from the one hypothesis or the other we shall arrive at widely 
differing results, which may be regarded, however, as the limits between 
which the facts lie. 

64. Assuming, in accordance with the first hypothesis, that the forma, 
tion of the continent in question involved a transfer of the earth's 
center of gravity towards the center of the continent, it will be of in- 
terest to compute the position of the disturbed sea surface with respect 
to spherical surfaces of the same radius and concentric about the origi- 
nal and disturbed centers of gravity, respectively. For this purpose 
we may use formulas (64) and (79). These are based on the assumption 
that the attracting mass is of uniform thickness. A more reasonable 
assumption is, perhaps, that the continents slope up rather rapidly, but 
not abruptly, from the sea shore, attaining their average height at no 

(163) 
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great distance inland; and more accnrate knowledge than we now 
possess might render it desirable to use some form of the more general 
equation (83), which takes account of variations in thickness of the 
attracting mass. As we can not expect, however, at present, to repre- 
sent the actnal shape of a continent very closely, and as we shall not 
attempt to estimate the attraction of the rearranged free water on itself, 
it will be best to confine attention to masses of uniform thickness, which 
give effects greater in general than the probable actual effects. 

65. For the angular radial extent and the relation of thickness and 
densities for the continent of Europe and Asia we take the following 
data substantially as they are given by Helmert in his Greodasie^, Part 
II, pp. 313, 314, viz : 

/?=arc of 380, 

3Ap/p^=4,000 meters 

=13,124 feet. 

The last expression is arrived at by taking for the average depth of 
the sea 3,438 meters, and for the average elevation of the continent 440 
meters, the density of the continental mass being assumed to be half 
the earth's mean density, i>r 2.8, and that of sea water 1. Thus we 
have a mass 3,438 meters thick, whose effective density is 2.8—1=1.8, 
and an additional mass 440 meters thick of density 2.8. These two are 
equivalent, so far as their potential to terms of the order we neglect is 
concerned, to a single mass of density 1.8 (or ^p^) and 4,120 meters 
thickness. Hence, in round numbers, the relation above. 

The value for the angular radius /3 of the continent is equivalent to 
about 2,600 miles, measured along the surface of the earth. It is the 
- * same radius assumed in Article XIY for the ice mass. This value is 
also very nearly that angular extent which a continent of uniform 
thickness must have to produce the maximum upheaval of water 
along its border (see section 23). 

The position of the disturbed surface relative to the two spherical 

surfaces concentric with the original and disturbed centers of gravity, 

;.• respectively, will be given at intervals of 10^ from the center of 

>i the continent to the point opposite, or 180^ from that center. The 

V' requisite values of I /n sin ^ for this jmrpose have been computed 

i from formulas (70) and (71), using the logarithms of ^i, ^2, ©tc, ^1, Jfc2, 

^'■' etc., given on page 69. ^o determine the slope of the disturbed surface 

\ . with respect to the spherical surfaces of reference (or the deflections of 

\ the plumb line) at intervals of 10^ from the center of the continent, 

pj • use has been made of equations (73) and (80), the differential' coef&cients 

'\ (iJ /da being computed from (74) and (76). Omitting the details of 

iV' the computation, the nature of which is readily apparent from the 

:j- ■ 

i' . * For lull title see page 86. 

(104) 
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equations referred to, the results are embodied in the follcwing table. 
The first column of the table gives angular distances increasing by in- 
crements of 10^ along a great circle from the center of the continent. 
The second and third columns giveforeach angular distance the elevation 
or depression of the disturbed surface relatively to the spherical sur- 
faces concentric with the undisturbed and disturbed centers of gravity 
of the earth, respectively, elevations being indicated b3' the plus sign 
and depressions by the minus sign. The fourth and fifth columns give 
for each angular distance the deflections of the plumb line or inclina- 
tions of the disturbed surface to the spherical surfaces of reference. 
The signs of these deflections are minus or plus according as the angle 
between the plumb line or normal at any point of the disturbed surface 
and the axis of the continent is greater or less than the angle between 
the radius vector of the satne point and the axis of the continent (see 
Fig. 7). 

2'able showing the superior limiting effects in disturbing the sea level attributable to the 

continent of Europe and Asia. 



Angular 

distance from 

center of 

continent. 


Elevation or depression of dis- 
turbed surface with respect to — 


Inclination of disturbed sur- 
face with reference to — 


Spherical sur- 
face concentric 
with undis- 
turbed center 
of gravity 
of earth. 


Spherical sur- 
face concentric 

with diS' 

tnrbed center 

of gravity 

of earth. 


Spherical sur- 
face concentric 

with undis- 
turbed center 

of gravity 
of earth. 


Spherical sur- 
face concentric 

with dis- 
turbed center 
of gravity 

of earth. 





Feet 


Feet. 


// 


// 


00 


+2,881 


+ 1,637 


- 0.0 


— 0.0 


30 


+2,812 


+1,587 


- 7.7 


- 5.6 


20 


* +2,588 


+1,419 


-17.7 


-13.5 


30 


+2, 148 


+1, 072 


-33.9 


-27.7 


40 


+1,099 


+ 146 


-61.3 


-43.5 


50 


^-h 444 


- 356 


-24.2 


-14.8 


60 


+ 111 


- 511 


-14.7 


- 4.1 


70 


- 104 


- 529 


-10.0 


+ 1.6 


80 


- 255 


- 471 


- 7.2 


+ 4.9 


90 


- 365 


- 365 


- 5.4 


+ 6.9 


100 


- 448 


- 232 


- 4.1 


+ 8.0 


no 


- 512 


- 87 


- 3.2 


+ 8.4 


120 


— 561 


+ 61 


- 2.4 


+ 8.2 


130 


- 600 


+ 200 


- 1.9 


+ 7.5 


140 


- 629 


+ 324 


- 1.4 


+ 6.5 


150 


— 650 


+ 426 


- 1.0 


. + 5.1 


160 


- 665 


+ 504 


- 0.7 


+ 3.5 


170 


— 673 


+ 552 


- 0.3 


+ 1.8 


180 


- 676 


+ 568 


- 0.0 


+ 0.0 


:^nll ±A 1 


I 
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66. The relations of the distnrbefl surface to the spherical sarAices 
of reference are shown in Fig. 7. The fall line circles in this diagram 
represent great circles of the spheres of reference on a scale of 2 50 oo ottoo • 
The position of the continent is indicated by the arc AB of the circle, 
whose center C is the nndistnrbed center of gravity of the earth. The 
disturbed center of gravity, or the center of the ojher circle, is at C 
The position of the disturbed surface is indicated by the dotted line. 
The radial distances of this line from the circles of reference and the 
distance between the centers G and C are exaggerated 2,500 times. 
QPR represents a normal to the disturbed surface at P. 

67. For the elevation of the disturbed surface at the immediate border' 
of the continent we find from the second of equations (07) and from 
(79) 1,380 feet or 400 feet, according as the sphere of reference is con- 
centric with the undisturbed or with the disturbed center of gravity of 
the earth. 

68. It should be remarked in this connection that the results given 
in the preceding table, and, indeed, all those in this section pertaining 
to the continent of Europe and Asia, are such as would exist under the 
assumed conditions if there was no other continent. The complete 
l)roblem requires the determination of the resultant action of all the 
continents at any point of the sea surface. For a method of determin- 
ing this resultant, when the components due to the several continents 
are known, the reader may be referred to Helmert's Geodasie. 

(1G6) 
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69. To determine the deflection of the plumb line at the level of the 
sea along the border of the continent, we observe that according to 
our assumed data the portion of the continental mass lying below the 
sea level is about 11,000 feet thick and has an effective density of 1.8, 
while the portion above the sea level is about 1,400 feet thick and has 
an effective density of 2.8. We may compute the deflection of the 
plumb line due to each of these portions by means of equation (77) and 
add the results to get the total deflection. For the lower mass we have 
for use in (77) 

/t= 11,000 feet, p=1.8, ^=11,000 feet; 

and for the upper mass 

^=1,400 feet, p=2.8, i?=0 feet. 

Hence the deflections due to the upper. and lower masses are 138'^ and 
34'', respectively, and their sum is 172''. This deflection is relative to 
a radius drawn to the undisturbed center of gravity of the earth. 
The deflection relative to a radius drawn to the disturbed center of 
gravity is by equation (80) 8" less, or 164". 

70. Under the conditions of the second hypothesis, which supposes 
the several radial element prisms of the earth's crust in a state border- 
ing on hydrostatic equilibrium, it is evident' that the disturbances of 
the sea level attributable to a continental mass must be of a low order. 
They must, in fact, be confined to terms of no higher order than those 
which have been neglected in our equations defining the position of the 
disturbed sea surface. The precise evaluation of these terms would offer 
difficulties practically insuperable in all cases, except those which present 
the simplest arrangement of densities in the element prisms of the earth's 
crust. We* may form a sufficiently definite idea of their smallness, 
however, byjeonsideriug the ideal question of the effect on the sea level 
of the radial transfer of a stratum of the earth's crust from some posi- 
tion below to some position above the sea surface. 

Let the stratum considered be a portion of a spherical shell ; let its 
border be circular and of angular radius /?. The effect on the sea level 
will obviously be greatest at the axis of the stratum. Hence, we only 
need to derive the change in potential at the sea level and at this axis 
due to the transferred stratum. Suppose the stratum raised to a height 
hi above sea level. Let its uniform thickness in this position be fi and 
its uniform density a« Then, r being the radius vector of any element 
of the mass and r' the radius vector of the sea surface, the potential 
of the stratum at its axis and at sea level will be 



sin eae 



2pi7r I f^dr / 

t/r'+A, J^ \/r2+r'2-2rr'cos^ 

(1G7) 
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This to terms of the first order inclusive is 

If we pappose the stratum in its original position is below the sea 
surface a distance i^, and that its uniform thickness and density in this 
position are tz and fhj respectively, the resulting potential at the same 
point of the sea surface is 

Since a^i=A^29 i* ^'y since the mass of the stratum in the two i>osi- 
tions is the same, thoagh its thickness and density may vary, the differ- 
ence of (97) and (98) is 



ir'pitiTt sin fi / Sfi+3h+ehi+Gh2 t,^t,-^2hi-^2h 






(99) 



This expression shows that a stratum of the earth's crust might be 
transferred through a considerable distance radially without materially 
affecting the sea level. 

Thus, for example, suppose 

^i=7i2=10,000 feet, ^,=^2=5,000 feet. 

Then the fraction in the parenthesis of (99) becomes -^. In other 
words, the potential at the sea surface under the center of the stratum 
would be only ^^ as great as that due to an uncompensated external 
stratum or internal vacuity of equal effective mass. 
In the case of the continent of Europe and Asia we have (section 69) 

7/1=0, fi=l,400 feet, pi=2.8. 

These in (97) give 

Ar'n sin § 1400 x 2.8/ 1+_L- 1— _ V («) 

^ I ^0000 60000 siu^y 

Likewise for use in (9S) we have 

A2=0, ^2=11,000 feet, p2=1.8. 

These give 

(16S) 
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The sum of (a) and {b) reduced to unit density is 



4r';r sin ^- 23 



''\ 3««» 9000 8lu§)" <") 



This is the potential to terms of the first order, inclusive of the whole 
continent at its center and at sea level, on the supposition that there 
is no defect in density of the underlying strata. But if we suppose a 
defect in density of 0.237 uniformly distributed through a depth of 
100.000 feet below the sea level, we must diminish the result (c) by 



The remainder is 

4r 



4r'^ P- 23720 /l — 1 ^ -\. 

2 V ''' 840 sin I J 

V /^- 23720 / — H ^^—77^, 

2 (,^«»%00sin§/ 



or about g-J^ of the potential due to the uncompensated mass. Under 
these conditions of compensation, therefore, the elevation of the sea at 
the center of the continent of Europe and Asia would be onl^^ about 10 
feet, whereas the first hypothesis would require an elevation of about 
2,900 feet (see section 65.) 

71. Notwithstanding the feeble effect a continental mass, whose radial 
elements are in a state approximating to hydrostatic equilibrium, would 
have in elevating or depressing the sea surface, it is conceivable that a 
considerable deflection of the plumb line might be produced by such a 
mass along its border. If, for example, we suppose the visible mass of 
the continent of Europe and Asia to be compensated by a defect in 
density of 0.18 uniformly dis1;ributed through a depth of 135,000 feet 
below sea level, the deflection of the plumb line along the border of this 
continent would be diminished by 117", leaving still a deflection (section 
69) of 172''— 117''=55" relative to the undisturbed sphere of reference. 
This deflection, 55", it should be observed, is the maximum possible 
value under the assumed conditions. If the degree of compensation as- 
sumed actually exists, it is probable that the real maximum deflection 
of the plumb line is much less than 65", since our calculation premises 
a vertical coast wall for the continent, whereas it presents for the most 
part only a moderately steep slope along the sea shore. 

XVII. LIST OF AUTHORS CONSULTED. 

72. The authors whose works have been specially consulted in the 
preparation of the preceding pages, and the full titles of their works, 
dates of publication, etc., are named below : 

Brans, Dr. Heinrich. Die Figur der £rde. Publication des konigl. preussischen 
geodatischen Institutes. Berlin, 1878. 
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- Clarke, Col. A. R. Geodesy. Oxford, 1880. 

Fischer, Dr. Pbilipp. Untersuchuageniiberdie Gestaltder Erde. Darmstadt, 1868, 
. Helmert, Dr. F. R. Die mathematischen und pbysikallsclien Tbeorieen der buberei> 
Geodasie. II Teil. Leipzig, 1884. 

Laplace. M^canique Celeste. Tome 3. Paris, 1878. 

Pratt, Jobn H. A Treatise on Attractions, Laplace's Functions, and tbe Figure of 
tbe Eartb. Fourth edition. London and New York, 1871. 

Stokes, G. G. Matbematical and Physical Papers, Vol. II, especially the paper On 
the Variation of Gravity at the Surface of the Earth. Cambridge University Pressy 
IS&X 

Thomson and Tait. Treatise on Natural Philosophy. Vol. I, Part II. Cambridge^ 
J1883. 

(170) 



' 



INDEX. 



ar Lake 

latics, mathematical feat- 

i in 

$ory, to relative position 

aces in a lake basin 

itinental glaciers 

ntinental masses 



Page. 
67 



mass M.. 

free water, 
ike basin.... 
glaciers 



13 

58 
60 
79 

19 
35 
58 
60 
79 
85 
17 

51 

of reference 19,51 

itthematical treatment of 

lems 13 

•iftless Area cited 17 

leroid of 19 



d, list of. 

surfaces in I^ke 

of earth and disturbing 



heroid 

3ference 

iluation of. 



19 
19 
37 



1 of glaciation 60 

•mergence 71 

fifect of ice mass 75 

>tion of form and position 

faces in lake basin 59 

>ntinental glaciers 61,79 

mtinent Europe and Asia 80 

nb-line, formulas for 47 

ike basin 60 

je mass 68 

ition of glaciation 70 

ent of Europe and Asia... 83, 85 
ea-level or geoid (see dis- 
use), von Drygalski on 78 

i in lake basin 58 

St 58,83 

II 59,61 

61 

;ase of continent Europe 

80 

rth's spheroid 19 

B masses 61 et seq. 

mass of Europe and Asia.. 80 
of gravity of earth 51 



Page. 

Disturbed surface, defined 19 

equations of. 40 

slope of. 47 

Driflless area, Chamberlin on the 13,17 

Drygalski, Erich von, memoir on Geoid- 

deformationen der Eiszeit cited... 73 

Effect of re-arranged free water 35 

of removal of water in lake basin 53 

of continental glaciers 60 

of continents on sea-level 79 

Elevation of disturbed surface 20,41,42 

Ellipsoid of reference, Clarke's elements of.. 19 

Equations of disturbed surface 40 

for disturbing mass of uniform thick- 
ness ;. 40 

for disturbing mass of variable thick- 
ness 53 

Equipotential surfaces, defined 12 

in lake basin 53 

Europe and Asia, continent of. 79-S5 

Evaluation of potential of disturbing 

mass 21 

of constants Fband Uo 37 

of definite integrals 7|, It 43 

Fischer, Philipp, treatise of 86 

Form, general of sea-level 15 

of disturbed surface 41 

of level surfaces in lake basin 58,59 

of sea-level as disturbed by assumed 

continental glaciers 66,70 

of sea-level as disturbed by a conti- 
nental mass 81,82 

Free water, attraction of 33 

Geodesy and geology as related to form 

and posilion of sea-level 15 

Geoid, position of, relative to spheroid 15 

deformation of, von Drygalski on 78 

Gilbert, G.K., request of, for investigation.. 13 

problem proposed by 17 

on slopes of lake beaches in New York.. 67 
Glaciers, continental, influence on sea- 
level 60 

assumed shape of 61 

slope near border of 62 

assumed volume of 64 

assumed mass of 65 

deflection of plumb-line on border of.... 68 
Glacial submergence 

discussion of, by Croll 71 

discussion of, by Pratt 71 

discussion of, by Heath 74 

discussion of, by Thomson 76 



87 



(171) 



88 



INDEX. 



Page. 
Glacial submcrgeuce— Contiuueil. 

discussion of, by von Dryjfalski 78 

table of data used by different authors 

in discussion of. 79 

Graphical illustration, 

of level surfaces in lake basin 59 

of shapes of assumed ice masses 63 

of disturbed surface due to alternation 

ofglaciation 70 

of disturbed surface due to continental 

masses 82 

Greenland, thickness of glaciers in 63 

Heath, D. D., on glacial submergence 74 

criticism of Croll and Pratt 75 

Helne^E., treatise cited 33 

Helmert, F. R., treatise cited 50,80,82 

Hypothesis relative to form of continental 

glaciers 61 

relative to nature of earth's crust 79 

Ice sheets, assumed slopes of surface and 

thickness of 61 

assumed volume of 64 

Illustrations, graphical, of level surfaces in 

lake basin 59 

of shapes of assumed ice masses 63 

of disturbed sea surface due to conti- 
nental glaciers 70 

of disturbed sea surface due to conti- 
nents 82 

Key to symbols used in this paper 9 

Lake Bonneville, beaches of. 17 

Lake Agifssiz, beaches of 67 

Lake basin, position of level surfaces in... 16.58 

L&ke Ontario, beach slopes near 67 

Laplace's functions 30 

Laplace, demonstration cited....*. 36 

Level surfaces, defined 18 

in lake basin 58 

slope of, in lake basin GO 

List of authors consulted 85 

Mass, of earth 19 

disturbing, defined 20 

in assumed ice-caps 64 

Mathematical statement of problem of this 

paper 18 

Mathematical symbols used in this paper, 

key to 9 

McGee, W J, on maximum synchronous 

glaciation 63 



Page. 
Nordenskjold on slope at border of Green- 
land ice-fields 63 

Ontario, slopes of beaches near Lake 67 

Plumb-line deflections, formulas for 47,51,52 

at border of assumed lake basin 60 

at border of assumed ice mass 68 

Potential of disturbing mass 19 

expressions for 24,33,35 

special values of 25 

degree of approximation of 26 

Position, relative, of disturbed surface 19,41 

of level surfaces in lake basin.s '. 58 

Pratt, Archdeacon J. H., on glacial sub- 
mergence 71-74,79 

Price, Bartholemew, treatise of, cited 27 

Slope of disturbed surface, formulas for...47, 51, 52 
Slope of upper surface of assumed ice 

masses 61,63 

Slope of glacial lake beaches 67 

Sphere of reference 19 

Spheroid of reference 15 

Clarke's elements of. 19 

Table of values showing 

relative position of level surfaces in 

lake basin 59 

slopes of assumed ice masses 62 

thickness of spherical shells of equal 
volume with assumed ice masses and 

equivalent lowering of sea-level 64 

position of disturbed surface in case of 

assumed ice masses 66 

minimum thicknesses of ice masses pro- 
ducing steep border slopes of dis- 
turbed suifice 68 

disturbance of sea-level due to assumed 

ice-cap 70 

superior limiting effects on sea-level at- 
tributable to continent of Europe 

and Asia 81 

Thomson, Sir W., on glacial submer- 
gence 71,76,79 

Undisturbed surface defined 19 

Upham, Professor Warren, on beaches and 

terraces of Lake Agassiz 67 

Variations in sea-level attributable to con- 
tinental glaciers 60 

attributable to alternation ofglaciation.. 69 

attributable to continental masses 79 

Works referred to 86 



(172) 



84 FORM AND POSITION OP THE SEA LEVEL. [bull. 48. 

This to terms of the first order inclusive is 

If we suppose the stratum in its original position is below the sea 
sarface a distance ^, and that its uniform thickness and density in this 
position are tz and p2, respectively, the resulting potential at the same 
point of the sea surface is 

Y 4r^8in^y 

Since piti^fhhy i» ©., since the mass of the stratum in the two posi- 
tions is the same, though its thickness and density may vary, the diflfer- 
ence of (97) and (98) is 

4r'M^ sin g/3^i+3f.+6/..+6/.._f.--f.+2/.,--2^.Y (99) 

V^ ^ 4r'sin^ J 



This expression shows that a stratum of the earth's crust might be 
transferred through a considerable distance radially without materially 
afifecting the sea level. 

Thus, for example, suppose 

^1=7*2=10,000 feet, ^,=^2=5,000 feet. 

Then the fraction in the parenthesis of (99) becomes 3^0 • I^ other 
words, the potential at the sea surface under the center of the stratum 
would be only ^^r ^s great as that due to an uncompensated external 
stratum or internal vacuity of equal effective mass. 
In the case of the continent of Europe and Asia we have (section 69) 

/«i=0, fi=l,400 feet, pi=2.8. 

These in (97) give 

4r';r sin § 1400x2.8/l+-L- 1 V («) 

2 ^ 20000 60000 Sin |y • 

Likewise for use in (9S) we have 

it 

*2=0, f2=l 1,000 feet, p2=1.8. 

These give 

ir'. Sin f llOOOx l.s/l-^-l-^-^^Y (,) 

' \ . i640sin^/ 

(16S) 



